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Abstract 

The  basic  concepts  for  implementing  today’s  Noise  Technology  Radar  (NTR) 
were  first  introduced  circa  1950.  However,  practical  implementability  was  severely 
limited  early  on  due  to  the  computational  intensity  of  analog  signal  correlation. 
Subsequent  advances  in  Digital  Signal  Processing  (DSP)  have  greatly  reduced  the 
computational  complexity  of  signal  correlation  using  a  Fast  Fourier  Transform  (FFT) 
and  increased  algorithm  efficiency  for  performance  optimization.  This  has  revived 
NTR  development  to  the  point  where  it  is  now  a  near-worldwide  activity,  with 
deployment  likely  to  occur  in  the  not  too  distant  future.  This  has  gotten  the  attention 
of  Radio  Frequency  Intelligence  (RFINT)  researchers  who  are  tasked  with  ensuring 
the  arsenal  of  Radio  Frequency  (RF)  situational  awareness  tools  remains  up  to  date; 
as  new  NTR  technology  emerges  so  does  the  need  for  non-cooperative,  non-matched 
filter  detection  and  exploitation. 

There  are  two  distinct  veins  of  ongoing  research  at  AFIT  involving  NTR  tech¬ 
nology,  including:  1)  NTR  developmental  work  aimed  at  improving  NTR  perfor¬ 
mance  from  a  cooperative  radar  detection  perspective  and  2)  NTR  exploitation  work 
aimed  at  assessing  NTR  exploitability  using  non-cooperative  detection  methods-the 
focus  of  this  research.  Activity  in  these  areas  is  based  on  AFIT’s  Noise  Network 
(NoNET)  radar  system  that  is  an  integrated  group  of  NTR  nodes  with  centralized 
processing  for  data  acquisition,  digital  processing,  and  multi-sensor  image  fusion. 
The  NTR  nodes  operate  in  Continuous  Wave  (GW)  mode  and  transmit  low  power 
signals  with  random  temporal-spectral  structure.  Thus,  the  signals  inherently  pos¬ 
sess  characteristics  that  are  commonly  associated  with  Low  Probability  of  Detec¬ 
tion/Intercept  (LPD/LPI)  waveforms. 
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As  adopted  for  this  research,  Quadrature  Mirror  Filter  Banks  (QMFB)  have 
been  effectively  used  for  non-cooperatively  detecting  structured  LPD/LPI  waveforms 
and  extracting  various  waveform  characteristics  (modulation,  pnlse  width,  band¬ 
width,  etc.).  The  goal  here  was  to  determine  the  effectiveness  of  QMFB  process 
for  extracting  features  from  AFIT’s  unstructured  NTR  signals  which  contain  no  in¬ 
tentional  temporal  or  spectral  modulation.  As  a  first  step,  this  was  done  through 
graphical  visual  assessment  of  energy-based  graphics  generated  directly  from  various 
QMFB  outputs-denoted  here  as  non-transformed  QMFB.  As  with  previous  research, 
the  QMFB  process  was  effective  at  detecting  structured  signal  features.  However, 
the  QMFB  process  performance  with  unstructured  NTR  signals  proved  to  be  more 
challenging  and  motivated  the  need  to  consider  transformed  QMFB. 

The  initial  approach  for  transformed  QMFB  was  based  on  AFIT  research 
by  Gronholz  and  Mims  (G-M),  who  successfully  demonstrated  a  method  for  non- 
cooperatively  detecting  LPD/LPI  LUtra  Wideband  (UWB)  communication  signals. 
The  G-M  process  uses  a  multi-channel  receiver  containing  a  bank  of  Narrowband 
(NB)  contiguous  Liters  that  span  the  UWB  signal  bandwidth.  The  filter  outputs 
are  synchronously  sampled  and  the  resulting  data  are  used  to  represent  the  signal’s 
Time-Frequency  (T-F)  response.  Given  common  LPD/LPI  characteristics  of  LTWB 
and  NTR  signals,  the  resultant  G-M  transformed  QMFB  with  NTR  signals  proved 
to  be  as  effective  as  what  was  demonstrated  with  LTWB  signals:  1)  For  collected  un¬ 
structured  NTR  signals,  there  was  minimal  discernible  difference  between  NTRqff 
and  NT Rqn  conditions,  and  detection  was  very  challenging,  and  2)  For  synthe¬ 
sized  structured  NTR  signals,  there  were  distinct  differences  between  NTRoff  and 
NT  Rqn  conditions  with  detection  occurring  through  visually-observed  features  that 
corresponded  to  the  induced  phase  features.  Overall,  the  research  goal  was  met- 
the  resultant  G-M  transformed  QMFB  is  useful  for  both  NTR  development  (desire 
to  maximize  LPD/LPI  potential)  and  exploitation  (desire  to  minimize  LPD/LPI 
potential)  research. 
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Characterization  of  Noise  Technology  Radar  (NTR) 
Signal  Detectability  Using  a  Non-Cooperative  Receiver 


I.  Introduction 

This  chapter  introduces  the  research  motivation,  objective  and  goal,  and  document 
organization.  The  research  motivation  is  provided  in  two  subsections,  with  one  pro¬ 
viding  the  Operational  Motivation  and  the  other  providing  Technical  Motivation. 
The  operational  motivation  provides  the  big  picture  concept  and  the  technical  moti¬ 
vation  identifies  key  concepts  that  served  as  the  research  foundation.  This  foundation 
was  used  to  formalize  the  research  objective  and  achieve  the  goal  of  advancing  the  ca¬ 
pability  to  non-cooperatively  detect  and  characterize  Noise  Technology  Radar  (NTR) 
signals.  This  chapter  concludes  with  an  overview  of  the  document’s  organization. 

1.1  Research  Motivation 

1.1.1  Operational  Motivation.  The  concept  of  noise  radar  was  first  intro¬ 
duced  circa  1950  [12, 18]  with  its  initial  implementability  primarily  limited  by  the 
computational  intensity  required  for  correlation  using  analog  circuits  [13].  Subse¬ 
quent  advances  in  computers,  technology,  and  software  have  revived  NTR  develop¬ 
ment  to  the  point  where  NTR  research  is  proliferating  throughout  the  world.  In 
particular,  advances  in  Digital  Signal  Processing  (DSP)  have  greatly  reduced  the  re¬ 
quired  computational  complexity  for  correlating  signals  using  a  discrete  Fast  Fourier 
Transform  (FFT).  At  the  same  time,  the  DSP  improvements  enable  implementation 
of  more  computationally-efficient  algorithms  and  performance  optimization.  There 
are  two  distinct  veins  of  ongoing  NTR  research,  including:  1)  NTR  developmen¬ 
tal  work  aimed  at  improving  NTR  performance  from  a  cooperative  radar  detection 
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perspective  [3-5, 10, 13, 17, 18,  25,  30-32]  and  2)  NTR  exploitation  work  aimed  at 
assessing  NTR  exploitability  using  non-cooperative,  “non-matched  filter”  detection 
methods  [11,27]. 

NTR  research  is  currently  underway  in  various  countries,  including:  the  United 
States,  Ukraine,  Poland,  United  Kingdom,  Canada,  Austria,  Montenegro,  France, 
Germany,  Italy,  Russia,  Sweden,  Norway,  China,  South  Korea,  and  India  [5, 10, 
13,17,18,25,30-32],  A  couple  specific  examples  of  global  research  abroad  include: 
1)  Swedish  researchers  seeking  to  improve  NTR  system  performance  using  beam¬ 
forming,  waveform  coding  for  Doppler  enhancement,  and  Synthetic  Aperture  Radar 
(SAR)  [2-4]  imaging;  and  2)  Polish  researchers  working  to  improving  NTR  system 
performance  using  impulse-based  sources  [14].  In  essence,  NTR  research  is  being 
conducted  worldwide  and  spans  a  vast  range  of  applications  from  SAR  imaging  to 
collision  avoidance  [2,18,19], 

The  near-worldwide  interest  in  NTR  system  development  and  its  likely  deploy¬ 
ment  in  the  not-too-distant  future  have  gotten  the  attention  of  researchers  in  the 
Radio  Frequency  (RF)  Intelligence  (RFINT)  community.  As  with  all  new  technology 
deployments,  the  RFINT  community  is  always  looking  forward  to  ensure  that  their 
arsenal  of  RF  awareness  tools  remains  up-to-date.  This  includes  the  capability  to 
non-cooperatively  exploit  identify  information  from  emerging  signals  of  interest. 

1.1.2  Technical  Motivation.  NTR  research  remains  an  ongoing  area  of  in¬ 
terest  at  the  Air  Force  Institute  of  Technology  (AFIT).  Specifically,  AFIT’s  NTR 
technology  is  being  developed  to  support  the  Noise  Network  (NoNET).  AFIT’s 
NoNET  radar  system  is  an  integrated  group  of  NTR  nodes  with  centralized  pro¬ 
cessing  for  data  acquisition,  digital  processing,  and  multi-sensor  image  fusion.  Of 
the  two  research  veins  noted  earlier,  NTR  developmental  work  has  been  the  primary 
focus  at  AFIT  [16,25,28,30].  This  is  due  to  the  low  power,  random  structure  of  the 
NTR  signal  which  result  in  electromagnetic  properties  that  cause  minimum  RF  inter- 
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ference  and  allow  the  NTR  system  to  share  frequency  bands  with  other  devices  [13]. 
Furthermore,  the  Continuous  Wave  (CW)  NTR  implementation  being  used  by  AFIT 
provides  range  resolution  that  is  on  par  with  conventional  CW  radar. 

The  emerging  NTR  signals  inherently  fall  into  the  class  of  Low  Probability  of 
Intercept  (LPI)  signals  given  that  they  are  relatively  low  power,  wide  bandwidth, 
and  have  frequency  variability  that  make  non-cooperative  detection  very  challeng¬ 
ing  [11,27].  Under  less  stricter  definitions,  they  are  also  Low  Probability  of  Detection 
(LPD)  signals  given  that  their  “presence”  is  difficult  to  determine  using  conventional 
power-based  detectors  such  as  radiometers  [8,21],  Thus,  the  LPD/LPI  characteris¬ 
tics  of  NTR  signals  being  developed  at  AFIT  provide  good  Electronic  Warfare  (EW) 
counter-countermeasure  capability  [13].  This  is  what  motivated  the  need  for  the 
second  vein  of  NTR  exploitation  research  as  conducted  under  this  research. 

The  LPD/LPI  research  in  [11,27]  provides  the  foundation  for  the  research  con¬ 
ducted.  Most  importantly,  these  previous  works  identified  the  Quadrature  Mirror 
Filter  Bank  (QMFB)  process  as  a  viable  approach  for  extracting  LPD/LPI  signal 
features.  The  QMFB  process  decomposes  the  LPD/LPI  signal  using  orthogonal 
wavelets  [6,11]  and  has  been  used  as  a  form  of  spread  spectrum  intercept  receiver  [6]. 
By  considering  the  QMFB  output  Time-Frequency  (T-F)  data  for  a  given  layer, 
graphic  representations  can  be  used  to  provide  relatively  good  estimates  of  param¬ 
eters  used  to  generate  structured  signals,  i.e.,  modulation  type,  bandwidth,  center 
frequency,  phase  modulation,  signal  duration,  etc.  [27].  It  is  important  to  note  that 
the  LPD/LPI  detection  work  in  [11,27]  was  based  on  structured  signals  containing 
distinct  modulation  features,  i.e.,  Frequency  Modulated  CW  (FMCW)  and  Binary 
Phase  Shift  Keyed  (BPSK)  signals.  Furthermore,  the  instantaneous  bandwidth  of 
the  non-cooperative  receiver  spanned  the  entire  bandwidth  of  the  FMCW  and  BPSK 
signals.  This  differs  considerably  from  the  research  conducted  here  which  is  based  on 
unstructured  LUtra  Wideband  (UWB)  NTR  signals  that  are  collected  using  a  Nar¬ 
rowband  (NB)  collection  receiver,  i.e.,  the  bandwidth  of  the  NTR  signal  (Wntr)  is 
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greater  than  the  collection  receiver  bandwidth  ( Wrf )  such  that  only  a  fraction  of  the 
NTR  signal’s  spectral  response  is  used  for  any  given  assessment  ( Wrf  ~  Wntr/W). 

Detection  and  estimation  of  unstructured  signals  such  as  those  used  for  NTR  is 
very  challenging.  Thus,  more  robust  detection  and  estimation  methods  were  sought 
under  this  research  using  QMFB  outputs  and  transformation  methods  investigated 
previously  by  Grohholz-Mims  (G-M)  for  UWB  communication  signals  [7,8,22,23]. 
The  G-M  process  is  a  non-cooperative  multi-channel  detection  method  that  was 
demonstrated  using  a  bank  of  NB  contiguously-spaced  filters  that  spanned  the  UWB 
communication  signal  bandwidth.  Given  the  characteristic  similarity  between  the 
UWB  communication  signal  and  the  UWB  NTR  signal  considered  under  this  re¬ 
search,  the  incorporation  of  the  G-M  process  was  a  reasonable  next  step  for  non- 
cooperative  NTR  signal  detection.  As  demonstrated,  the  G-M  process  exploits  in¬ 
herent  temporal  and  spectral  differences  using  four  detection  methods.  Each  of  these 
four  methods  were  assessed  for  the  NTR  signal  of  interest,  (SOI). 

1.2  Research  Objective  and  Goal 

The  overall  research  objective  was  to  collect,  detect,  and  assess  NTR  signal 
characteristics  using  qualitative  graphical  visualization.  This  objective  was  formal¬ 
ized  relative  to  the  two  veins  of  research  being  supported,  including  development  (de¬ 
sire  to  maximize  LPD/LPI  potential)  and  exploitation  (desire  to  minimize  LPD/LPI 
potential).  The  goal  was  to  improve  upon  existing  QMFB  signal  processing  and 
graphical  visualization  methods  to  enable  non-cooperative  detection  of  NTR  signals. 
The  resultant  analysis  graphics  are  similar  to  what  was  envisioned  for  an  input  to  a 
human  interface  device  to  enable  an  operator  to  detect  NTR  signals  (declare  NTR 
presence).  The  detection  approach  used  a  non-cooperative  collection  receiver  and 
empirical  pair-wise  assessment  of  graphics  obtained  for  NTRoff  and  NT  Ron  col¬ 
lection  conditions.  The  emergence  of  distinguishable  feature  characteristics  and/or 
differences  between  the  two  conditions  constituted  NTR  signal  detection. 
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1.3  Document  Organization 

This  document  includes  five  chapters  and  three  appendices.  This  section  out¬ 
lines  what  is  provided  in  each  chapter.  Chapter  II  provides  relevant  technical  back¬ 
ground  information  on  three  major  concepts  used  to  support  the  research,  including 
1)  AFIT’s  NTR-Based  NoNET,  2)  QMFB  processing,  and  3)  G-M  processing. 

Chapter  III  provides  the  overall  approach,  process,  and  techniques  used  to  visu¬ 
ally  assess  and  characterize  the  NTR  signal.  A  detailed  description  of  the  collection 
process  and  NTR  configurations  is  first  presented.  This  is  followed  by  descriptions 
of  various  post-collection  processes  that  were  implemented  in  MATLAB®. 

Chapter  IV  provides  qualitative  graphical  results  and  assessments  using  the 
process  outlined  in  Chapter  III.  Verification  of  the  QMFB  process  is  first  provided, 
followed  by  the  step-by-step  application  of  processes  outlined  in  Chapter  III  with 
key  down-selection  details  highlighted  where  appropriate.  Finally,  the  overall  results 
are  presented  for  selected  scenarios;  the  amount  of  available  graphical  data  is  quite 
immense  and  only  a  fraction  of  it  is  presented  to  ensure  that  a  succinct  textual 
presentation  is  made  on  NTR  detectability. 

Chapter  V  concludes  the  thesis  and  provides  an  overall  summary  of  the  research 
motivation,  goals,  techniques,  results,  and  follow-on  recommendations.  This  section 
is  followed  by  three  appendices  that  that  provide  supplementary  graphical  results 
and  the  MATLAB®  code  to  implement  select  post-collection  processes. 
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II.  Background 

This  chapter  presents  relevant  technical  background  information  on  three  major 
concepts  and  techniques  used  to  conduct  this  research,  including:  1)  AFIT’s  NTR- 
Based  Noise  Network  (NoNET),  2)  Quadrature  Mirror  Filter  Bank  (QMFB)  pro¬ 
cessing,  and  3)  Gronholz-Mims  (G-M)  processing.  As  described  in  Section  2.1,  the 
primary  Signal  of  Interest  (SOI)  is  the  Noise  Technology  Radar  (NTR)  signal  being 
developed  for  use  in  AFIT’s  NoNET.  The  non-cooperative,  “non-matched  filter”  de¬ 
tection,  estimation,  and  exploitation  of  NTR  signal  features  are  addressed  using  the 
QMFB  process  described  in  Section  2.2.  Improvement  in  NTR  detectability  through 
graphical  visual  assessment  is  addressed  using  the  G-M  process  in  Section  2.3.  For 
both  non-transformed  QMFB  and  G-M  transformed  QMFB  processes,  NTR  sig¬ 
nal  detectability  is  assessed  using  comparative  analysis  of  graphics  generated  under 
NTRoff  and  NT  Ron  conditions. 

2.1  AFIT’s  NTR-Based  Noise  Network  (NoNET) 

The  NTR  signals  considered  for  this  research  were  generated  using  one  of  the 
NTR  systems  (nodes)  in  AFIT’s  Noise  Network  (NoNET)  operating  a  monostatic 
configuration  in  Figure  2.1.  The  NoNET  was  first  established  at  AFIT  in  2009  [30] 
and  is  based  on  the  developmental  work  by  Dr.  Narayanan  at  Pennsylvania  State 
University  [24].  Since  that  time,  work  has  continued  at  AFIT  to  improve  NoNET 
monostatic  and  multistatic  radar  performance  [25].  Currently,  there  are  two  distinct 
veins  of  research  involving  NoNET  research,  including  1)  NTR  developmental  work 
to  improve  NoNET  performance  from  a  radar  perspective  [16, 28]  and  2)  NTR  ex¬ 
ploitation  work  to  assess  NoNET  exploitability-the  focus  of  this  research.  The  NTR 
development  component  aims  to  improve  cooperative  transmit-reflect-receive  signal 
detectability  while  the  exploitation  component  aims  to  improve  the  non-cooperative 
detectability  of  transmit-only  NTR  signals. 
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Figure  2.1:  Functional  diagram  of  AFIT’s  NTR-Based  Noise  Network 
(NoNET)  radar  system  [28] 

2.1.1  NoNET  NTR  Implementation.  As  detailed  in  Figure  2.2  and  used 
for  this  research,  each  NoNET  node  can  be  operated  as  an  independent  monostatic 
NTR.  The  basic  NTR  system  is  comprised  of  three  main  components,  including  a: 
1)  Transmitter,  2)  Receiver  and  3)  Receiver  Processor.  A  majority  of  the  NTR 
system  complexity  resides  in  the  receiver  processor,  which  is  where  correlation  is 
used  to  detect  target  responses. 


NTR  Receiver 


Figure  2.2:  Implementation  diagram  of  typical  monostatic  NTR  node 
in  AFIT’s  NoNET  radar  system. 


As  currently  configured,  the  NTR  operates  as  a  CW  radar  that  ideally  trans¬ 
mits  a  purely  random  signal  produced  by  a  solid-state  thermal  noise  source.  The 
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thermal  noise  response  can  be  characterized  as  bandlimited  white  Gaussian  noise 
given  that  its  Power  Spectral  Density  (PSD)  is  nearly  uniform  over  a  finite  band¬ 
width  [25,26].  The  noise  source  output  is  effectively  bandpass-filtered  using  a  se¬ 
ries  of  low-pass  and  high-pass  filters.  The  resultant  NTR  transmit  bandwidth  is 
W nt r  ~  400  MHz  [30].  The  filtered  signal  is  power-divided  to  create  the  transmit¬ 
ter  output  signal  ( Tx  Out),  Sr(t),  which  is  fed  to  a  high  gain  log-periodic  antenna 
for  over-the-air  transmission.  The  other  output  is  fed  to  the  receiver  processor  and 
serves  as  the  transmitter  reference  ( Tx  Ref),  SRef(t).  The  Tx  Out  signal  propagates 
through  the  channel,  reflects  off  a  target  (s),  and  is  received  by  another  high  gain 
log-periodic  antenna.  The  received  input  signal  ( Rx  In),  s^ft),  is  filtered  to  ex¬ 
tract  desired  frequency  components  and  passed  to  the  receiver  processor  for  digital 
correlation  with  a  stored  copy  of  SRef(t)  [25,30]. 

The  correlation  process  in  the  NTR  receiver  processor  effectively  distinguishes 
the  target  return  signal  (NTR  noise  signal)  from  nndesired  random  noise  responses 
induced  by  the  channel  and  present  in  the  received  signal.  Given  that  the  transmit¬ 
ted  signal  response  is  known,  SRef(t)  =  sr(t)  and  cross-correlation  can  be  used  to 
determine  the  target  presence  in  the  return  signal.  This  is  implemented  using 

CReffiij)  =  f  SRef(t)sR(t  +  r)dt,  (2.1) 

Jo 

where  T  is  the  time  duration  of  the  recorded  transmitted  signal  at  each  frequency  and 
t  is  the  time  delay.  The  function  in  (2.1)  provides  a  measure  of  “likeness”  between 
the  signals  being  correlated.  When  SRef(t )  is  present  in  SR(t),  the  system  produces 
a  high  cross-correlation  value  and  a  “thumbtack”  like  graphical  response  which  is 
referred  to  as  auto-correlation.  At  sufficiently  high  SNR,  this  peak  response  rises 
above  the  noise-only  response  and  the  presence  of  an  NTR  signal  is  detected. 

The  correlation  in  (2.1)  is  functionally  implemented  digitally  for  computational 
efficiency.  The  notations  SRef[n]  and  s r [n]  are  introduced  here  to  represent  sampled 


versions  of  SRef(t )  and  sR(t),  respectively.  Specifically,  the  digital  implementation 
includes  1)  converting  SRef[n }  and  ,s/?[n]  to  the  spectral  domain  using  a  discrete 
Fourier  Transform  (DFT),  2)  multiplying  the  resultant  Fourier  components,  and 
3)  taking  the  inverse  discrete  Fourier  Transform  (1DFT)  returns  the  resultant  to  the 
sampled  time  domain.  DFT  and  1DFT  are  computed  using  FFT.  This  is  analytically 
expressed  as 


CRef,R[n]  =  IDFT  {SRef(k)S*R(k)}  ,  (2.2) 

where  SRef[k]  and  Sr [k]  the  DFT  of  SRef[n]  and  srt[n],  and  *  is  the  complex  conju¬ 
gate.  With  a  known  transmitted  signal  and  correlation,  the  NTR  system  is  a  very 
robust  and  effective  process  to  detect  a  returned  target  signal. 

2.1.2  NTR  Signal  Detection.  There  are  various  terminologies  used  through¬ 
out  the  NTR  community  that  refer  to  specific  implementations  of  NTR  systems,  two 
of  which  include  Random  Signal  Radar  (RSR)  and  Random  Noise  Radar  (RNR)  [9, 
11,  15].  RNR  refers  to  a  system  using  a  non-modulated  unstructured  noise  sig¬ 
nal  [11,15]  while  RSR  refers  to  a  system  using  a  modulated  structured  random  noise 
signal  [9,11].  Collectively  considering  all  of  these  related  works,  NTR  developmental 
research  appears  to  be  a  global  activity.  When  considering  the  non-cooperative  ex¬ 
ploitation  of  NTR  the  work  appears  to  be  limited  to  a  few  researchers  [11,27].  These 
non-cooperative  activities  generally  classify  NTR  signals  as  being  Low  Probability  of 
Detection  (LPD)  and  Low  Probability  of  Intercept  (LPI)-complex  signals  that  are 
difficult  to  detect  and  which  possess  covert  characteristics. 

2.1.3  LPD/LPI  Signal  Detection.  The  work  in  [11,27]  considered  four 
different  detection  methods  for  LPD/LPI  signals,  including:  Wigner-Villc,  Choi- 
Williams,  QMFB,  and  cyclostationary  processes  [27].  Based  on  this  earlier  research, 
QMFB  process  emerged  as  the  technique  that  provided  the  most  potential  for  reveal- 
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ing  signal  features  in  “structured”  signals,  i.e.,  signals  containing  distinct  features 
in  a  given  signal  domain(s)  (time,  frequency,  etc.).  The  QMFB  process  decomposes 
the  LPD/LPI  signal  using  orthogonal  wavelets  [6,11]  and  has  been  used  as  a  form 
of  spread  spectrum  intercept  receiver  [6].  By  considering  the  QMFB  output  Time- 
Frequency  (T-F)  data  for  a  given  QMFB  layer,  graphic  visualization  can  be  used 
to  provide  relatively  good  estimates  of  the  structured  signal  parameters,  i.e.,  mod¬ 
ulation  type,  bandwidth,  center  frequency,  phase  coding,  and  signal  duration  [27]. 
The  previous  research  in  [11,27]  was  limited  to  considering  structured  signals,  i.e. 
FMCW  and  BPSK  signals. 

Laver  1  Layer  2  Layer  3 


Figure  2.3:  QMFB  process  overview  [11,27]. 


2.2  QMFB  Processing 

The  QMFB  process  sequentially  “breaks  down”  the  in-phase  and  quadrature- 
phase  (I&Q)  components  of  the  input  signal  using  modified  sine-shaped  filters  in 
a  series  of  Quadrature  Mirror  Filter  (QMF)  pairs  and  layers  as  illustrated  in  Fig¬ 
ure  2.3.  The  output  of  each  layer  is  T-F  data  for  each  QMF  pair,  with  the  amount 
of  QMF  pairs  increasing  by  a  power  of  two  with  each  increasing  layer  number.  As 
illustrated  in  the  bottom  of  Figure  2.3,  the  QMFB  layer  outputs  can  be  represented 
as  rectangular  T-F  “tiles”  with  the  tile  dimensions  (T-F  resolution)  varying  as  a 
function  of  layer  number.  For  this  research,  the  QMFB  process  output  for  layer  n  is 
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denoted  by  Qn  and  referred  to  as  the  non-transformed  QMFB  response.  Represen¬ 
tative  graphics  for  a  QMFB  processed  FMCW  waveform  are  provided  in  Figure  2.4 
for  illustration.  These  particular  graphics  are  based  on  Layer  ffQ  output  data  (Q6) 
and  were  generated  using  MATLAB®  graphic  functions  -contour,  mesh,  surface  and 
Pseudo-Color  (denoted  henceforth  as  P color). 
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(a)  Q 6  Contour  Graphic. 


(c)  Qq  Surface  Graphic. 


(b)  Q6  Mesh  Graphic. 


Time  (sec) 

(d)  Q 6  Pcolor  Graphic. 


Figure  2.4:  Four  graphical  representation  (as  indicated)  of  the  Non- 
Transformed  QMFB  (Q6)  for  Layer  ff 6  NT  Ron  condition  using  Nq  = 
213  samples. 


With  each  increase  in  layer  number,  there  is  a  corresponding  trade-off  in  T-F 
resolution  [27],  i.e. ,  frequency  resolution  increases  along  the  frequency  dimension  (F) 
as  the  signal  is  decomposed  (filtered)  at  each  stage  into  low  and  high  frequency  com¬ 
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ponents.  Correspondingly,  the  time  resolution  decreases  along  the  time  dimension 
(T)  given  the  inversely  proportional  T-F  relationship  [6].  The  time  resolution  (At) 
and  frequency  resolution  (A/)  are  given  by  [27] 


A/ 

At 


fs  =  h_ 

2  (2^-1  -  1)  Nf  ’ 

2  nl  2Nl 

{,  (2^-'  -  1)  =  JJh  ’ 
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where  Nl  is  the  total  number  of  QMFB  layers,  Np  is  the  number  of  tiles  along  the 
frequency  dimension,  Np  is  the  number  of  tiles  along  the  time  dimension,  and  fs  is 
the  sampling  frequency.  Consistent  with  the  T-F  inverse  proportional  relationship, 
lower  numbered  layers  have  higher  time  resolution  and  lower  frequency  resolution 
when  compared  to  higher  numbered  layers  which  have  reduced  time  resolution  and 
higher  frequency  resolution. 

The  frequency  decimation  in  the  QMFB  process  is  implemented  here  using  a 
“modified  sine”  Finite  Impulse  Response  (FIR)  filters  having  a  response  given  by  [6] 
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where  S  is  the  scaling  variable,  C  is  the  compression  variable,  N  is  the  number  of 
coefficients,  and  w  [n]  is  a  Hamming  window  weighting.  The  Hamming  window  is  used 
to  suppress  the  effects  of  Gibb’s  phenomena  resulting  from  truncation  of  n  [6].  These 
particular  filters  have  a  flat  pass  region  and  produce  the  maximum  amount  of  desired 
signal  energy  within  each  tile  while  rejecting  other  unwanted  responses.  Based  on 
research  in  [6],  the  ideal  filter  coefficients  for  N  =  512  are  C  =  1.99375872338059 
and  S  =  1.00618488680080  which  provide  a  set  of  “nearly  orthogonal  filters  with 
cross-correlation  of  less  than  .001”  [6,11,27]. 

There  are  some  special  conditions  to  consider  when  implementing  QMFB  pro¬ 
cess.  For  example,  the  number  of  input  samples  to  the  QMFB  process  (Nq)  must 
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satisfy  Nq  =  2n  for  n  a  positive  integer.  This  produces  Nl  QMFB  layers: 

Nl  =  log2  (Nq)  .  (2.6) 

For  a  collected  signal  sample  size  of  Ns  7^  2n,  the  collected  signal  can  either  be 
parsed  to  a  power-of-two  producing  Nl  layers  per  (2.7)  or  zero-padded  and  extended 
to  a  power-of-two  to  produce  NL  layers  per  (2.8). 

Nl  =  floor[log2(Ns )]  .  (2.7) 

Nl  =  ceil[log2(N s)]  ■  (2.8) 

The  resultant  matrices  (Qn)  from  QMFB  process  may  be  used  for  signal  de¬ 
tection  (establishing  presence)  and/or  assessment  of  signal  features  as  captured  in 
the  Qn  T-F  data.  Generally,  Qn  matrices  contain  different  information  that  can  be 
visualized  through  various  graphical  representations.  The  amount  of  information  re¬ 
vealed  is  a  function  of  the  layer  number,  n,  and  the  specific  graphic  used  to  visualize 
the  transformed  data  (surface,  contour,  etc.) 

2.3  Gronholz-Mims  (G-M)  Processing 

As  adopted  here,  the  Gronholz-Mims  (G-M)  process  was  developed  in  [7,  8, 
22, 23]  to  enable  Narrowband  (NB)  processing  and  exploitation  of  Ultra  Wideband 
(UWB)  Communication  Signals.  The  main  goal  in  the  G-M  process  is  to  exploit 
potential  temporal  and  spectral  differences  in  NB  signal  responses  by  transforming 
T-F  data  to  produce  four  alternate  data  matrices,  including:  Temporal- Temporal 
Matrix  (TTM),  Cross- Temporal  Matrix  (CTM),  Spectral-Spectral  Matrix  (SSM), 
and  the  Cross- Spectral  matrix  (GSM).  The  G-M  process  is  depicted  in  Figure  2.6 
which  shows  the  complex  Channelized  Data  Matrix  (CDM)  as  the  input. 


13 


For  previous  G-M  demonstrations  in  [7,8,22,23]  the  input  CDM  was  generated 
using  the  channelized  receiver  architecture  illustrated  in  Figure  2.5.  In  this  case, 
the  CDM  effectively  represents  the  MqxNq  T-F  response  of  the  input  signal,  where 
r(t)  is  the  UWB  communication  signal,  Mq  is  the  number  NB  receiver  channels, 
Nq  is  the  number  of  synchronous  time  samples  taken  for  each  channel,  and  Hrf  is 
the  input  UWB  communication  signal  bandwidth.  The  receiver  is  a  bank  of  equal 
bandwidth,  NB  bandpass  filters  with  different  center  frequencies  that  collectively 
span  the  full  bandwidth  of  the  UWB  input  signal.  For  general  application,  the  G-M 
process  in  Figure  2.6  can  be  applied  using  any  input  matrix  representing  T-F  data. 
The  process  flow  for  G-M  transformation  per  Figure  2.6  is  described  in  the  next 
paragraph. 


Figure  2.5:  G-M  UWB  Channelized  Receiver  [7,8,22,23]. 
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Figure  2.6:  Gronholz-Mims  (G-M)  process  overview  [7,8,22,23]. 
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In  order  to  generate  the  TTM  and  corresponding  CTM  matrices,  the  input 
CDM  is  first  transformed  using  an  A^FFT-point  IFFT  on  a  column-by-column  basis 
( column-wise )  to  produce  the  TTM.  If  NIFFT  >  Mq  zero-padding  can  be  used  and  if 
NiFFt  <  Mq ,  the  higher  frequency  elements  may  be  neglected.  The  corresponding 
column-wise  and  row-wise  CTM  is  easily  generated  using  [7, 8, 22, 23] 


CTMcoi  = 


N] 


IFFT 


TTM 


H 


TTM 


(2.9) 


or 


CT  MRow  — 


1 

Njppt 


TTM  ■  TTMh 


(2.10) 


where  the  superscripted  H  denotes  the  Hermitian,  or  conjugate  transpose  operation. 
Also,  the  CTM  can  be  created  using  the  rows  of  TTM. 

The  SSM  and  corresponding  CSM  matrices  are  generated  using  a  similar  pro¬ 
cess,  with  the  input  CDM  is  first  transformed  using  an  NFFF-pomt  FFT  on  a  row-by- 
row  basis  ( row-wise )  to  produce  the  SSM.  If  NFFF  >  Nq  zero-padding  may  be  used 
and  if  NFFF  <  Nq,  the  higher  time  elements  may  be  neglected.  The  corresponding 
column-wise  and  row-wise  CSM  is  easily  generated  using  [7,8,22,23] 


CTMCol  =  — 

iv  fft 


ssmh  ■  SSM 


(2.11) 


or 


CSMRow  —  — - 

iv  fft 


SSM  ■  SSMh 


(2.12) 


The  resultant  matrices  from  (2.9)  through  (2.12)  may  be  used  for  signal  de¬ 
tection  (establishing  presence)  and/or  assessment  of  signal  features  as  captured  in 
the  input  CDM  T-F  data.  Generally,  the  four  resultant  G-M  matrices  contain  differ¬ 
ent  information  that  can  be  visualized  through  graphical  representation  using  P color 
graphic  function  in  MATLAB®.  The  amount  of  information  revealed  is  a  function  of 
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where  the  information  exists  (temporal  domain,  spectral  domain,  or  cross-correlation 
thereof). 

Summary 

This  chapter  presented  the  relevant  background  information  on  key  technical 
concepts  required  to  conduct  this  research.  This  included  a  summary  of  1)  AFIT’s 
NTR-Based  Noise  Network  (NoNET),  2)  Quadrature  Mirror  Filter  Bank  (QMFB) 
Processing,  and  the  3)  Gronholz-Mims  (G-M)  Processing.  Sufficient  references  are 
provided  in  respective  sections  if  more  detailed  information  is  required. 
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III.  Methodology 

This  chapter  provides  the  research  methodology  used  for  assessing  NTR  signal  de¬ 
tectability  by  applying  the  process  illustrated  in  Figure  3.1.  This  process  was  used 
to  support  the  technical  contribution  of  this  research  in  two  aspects: 

1.  This  research  represents  the  first  known  assessment  of  NTR  detectability  us¬ 
ing  the  (Quadrature  Mirror  Filter  Bank)  QMFB  process  and  the  Previous 
graphics  that  are  commonly  used  [27].  Thus,  the  Previous  label  in  Figure  3.1 
simply  indicates  that  the  graphic  visualization  is  the  same  as  used  previously 
with  structured  signals.  The  importance  of  applying  the  end-to-end  process  in 
Figure  3.1  with  NTR  signals,  i.e.,  from  NTR  signal  collection  through  post- 
collection  processing  and  visualization  with  Previous  graphics  should  not  be 
overlooked. 

2.  The  dashed-box  with  the  Current  label  in  Figure  3.1  represents  implementa¬ 
tion  of  multi-channel  processing  and  graphic  visualization  used  previously  by 
Gronholz  [7,  8]  and  Mims  [22,  23]  to  qualitatively  assess  the  detectability  of 
Ultra  Wideband  (UWB)  communication  signals.  Of  direct  relevance  to  this 
research,  these  previous  methods  are  based  on  (Narrowband)  NB  assessment 
of  UWB  signals-the  approach  adopted  here  for  NTR  signals.  To  identify  the 
application  of  these  previous  methods  with  results  generated  here,  the  (G-M) 
notation  is  introduced  here  and  used  henceforth  in  the  document. 

This  chapter  includes  six  major  subsections:  1)  Pre-Conditioning  Operations, 
2)  RF  Signal  Collection,  3)  Simulated  FMCW  Signal,  4)  Experimental  NTR  Signal, 
5)  Post-Collection  Processing,  6)  Graphic  Visualization,  and  6)  Wealth  of  Available 
Data.  Section  3.2  provides  an  understanding  of  the  collection  system,  Radio  Fre¬ 
quency  System  Intercept  and  Collection  System  (RFSICS).  Next,  Section  3.3  and 
Section  3.4  describe  the  SOI  FMCW  and  NTR  with  detailed  descriptions  of  the 
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three  NTR  configurations.  Section  3.5  goes  into  the  main  research  area  covering  the 
post-collection  and  post-QMFB  processing,  minus  the  graphic  visualizations  which 
are  covered  in  the  next  section. 


Previous  Current 

Figure  3.1:  Overview  of  process  used  to  assessing  NTR  detectability. 


3.1  Pre-Conditioning  Operations 

This  section  provides  details  of  pre-conditioning  operations  that  were  applied  to 
various  complex  data  matrices  throughout  the  overall  process  detailed  in  Figure  3.1. 
The  operations  are  presented  in  no  particular  order  using  the  notation  Q(j,  k)  to 
represent  the  jth  row  and  kth  column  of  an  arbitrary  Complex  Matrix  having  Nj  rows 
and  Nk  columns.  As  ultimately  used  for  the  research  and  detailed  in  Section  3.6.1 
and  Section  3.6.2,  the  pre-conditioning  operations  can  be  applied  in  any  order  to: 

1.  Any  of  the  complex  matrices  shown  in  Figure  3.1,  including  the  direct  QMFB 
output  Qn  and  the  G-M  transformed  Q%M:  Qttmi  Qctmj  QssmiQcsm- 
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2.  Any  resultant  output  matrix  from  previous  pre-conditioning  operations,  i.e., 
the  operations  can  be  sequentially  applied.  For  example,  the  input  matrix 
Real  [Q]  may  be  first  row-wise  magnitude  normalized  using  (3.3),  followed  by 
row-wise  cross-correlated  using  (3.10),  and  converted  to  dB  using  (3.1)  before 
graphics  are  generated. 

The  available  pre-conditioning  operations  include: 

1.  Available  Inputs:  Q  (C),  Real[Q]  (M),  Imag[Q]  (I),  and  Q\  (Magnitude). 

2.  Converting  to  Energy  in  dB: 

E%  =  10Log10(\Q\2)  (3.1) 


3.  Global  Maximum  Magnitude  Normalization: 

QU,  k)  =  QU,  k)/max  [|  Q\\  Vj,  k  (3.2) 

4.  Row-Wise  and  Column-Wise  Maximum  Magnitude  Normalization: 


QU,  ■)  =  QU,  0 /max  [|  Q(j,  :)|]  Vj 

(3.3) 

Q(:,k )  =  Q(:,k)/max[\Q(:,k)\\  \/k 

(3.4) 

5.  Row-Wise  and  Column-Wise  Power  Normalization: 

pqH)  QU, k )2  ■, 

k= i 

QU,  k )  =  x  Q(j,  k )  ,  Vj,  k 

(3.5) 

Ni 

W«FEQ(J^)2,  VA: 

J  3= 1 

QU,  k )  =  x  QU,  k)  ,  Vj,  k 

(3.6) 

6.  Row-Wise  and  Column-Wise  Amplitude  Normalization: 

AqU)  =  max[\QU,-)\\  ,  Vj 

QU,  k)  =  QU,  k)/AQU)  ,  Vj,  k 

(3.7) 
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AQ{k)  =  max[\Q(:,k)\]  ,  VA: 

QU,k)=QV,k)/AQ(k)  ,  Vj,  k 

7.  Thresholding/Flooring: 

Calculate  :  tQ  =  min  [Q] 

Initialize  :  QF(j,k)  =  Q(j,k)  ,  Vj,  k  (3.9) 

Floor  :  QF(j,  k)  <  tQ  =>  QF(j,  k)  =  tQ  ,  Vj,  k 


8.  Thresholding/Flooring,  Pair-Wise  Associated  NTRoff  and  NTRon- 

Calculate  :  tq  =  max  {min  [  Qqff }  >  mi'ri  [Qoiv]} 

Initialize:  QoFFUik)  =  QOFF(j,k )  ,  Vj,  k 

Initialize  :  QoMk)  =  QONU,k)  ,  Vj,Ar  (3-10) 

Floor  :  QoFF(j,  k)  <  tQ  =*►  QoFF(j,  k)  =  tQ  ,  Vj,  k 

Floor  :  Q%N(j,  k)  <  tQ  =>  Qav(A  k)  =  tQ  ,  Vj,  k 

9.  Row-Wise  and  Column-Wise  Cross-Correlation  [7,8]: 

Qg”1  =  tJt  (Qh  ■  Q) 

where  superscripted  H  denotes  Hermitian  or  conjugate-transpose,  Aj-  is  the 
number  of  rows  and  IFFT  operations,  and  Aj  is  the  number  of  columns  and 
FFT  operations. 


3.2  RF  Signal  Collection 

The  RFSICS  is  used  to  collect  NTR  signals  using  the  procedures  detailed  in 
Appendix  A  of  [29].  A  description  of  key  RFSICS  characteristics  is  provided  in 
Section  3.2.1.  The  RFSICS  was  used  for  the  three  NTR  collection  configurations 
detailed  in  Section  3.4.1,  Section  3.4.2,  and  Section  3.4.3.  For  all  configurations 
the  RFSICS  was  first  set  up  to  collect  with  the  NTRon  condition  to  enable  set¬ 
ting  of  the  Analog-to-Digital  converter  (ADC)  dynamic  range  to  avoid  clipping  the 
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collected  signals.  The  NTR  system  was  then  switched  off  and  collected  signals  for 
the  NTRoff  condition.  Thus,  for  all  comparative  NT  Roff  /  NT  Ron  assessment, 
the  variability  due  to  collection  system  bias  is  minimal  given  that  the  RFSICS  set¬ 
tings  were  identical  for  the  pair  of  NTRoff  and  NT  Ron  signals  being  assessed. 
Comparing  the  collection  conditions  is  a  start  to  determine  if  the  NTR  systems  is 
detectable  and  reveal  any  information  that  is  inherent  to  the  transmitter.  Follow¬ 
ing  collection,  the  signal  was  converted  from  the  Agilent®  proprietary  “capture”  file 
(*.cap)  to  standard  MATLAB®  (*  .mat)  format  for  post-collection  processing  [29]. 
The  resultant  collected  data  sequence,  denoted  as  {/?(&)},  contains  (Ns)  collected 
complex  I&Q  samples. 

3.2.1  RFSICS  Characteristics.  The  RFSICS  is  an  Agilent®-based  re¬ 
ceiver  system  that  AFIT  uses  for  RF  signal  collection  and  experimentation  [1], 
The  RFSICS  has  a  Narrowband  (NB)  instantaneous  collection  bandwidth  (Wrf) 
of  Wrf  =  36  MHz.  This  is  coupled  with  Wideband  (WB)  search  capability  across 
the  range  of  fcE  [0.02,  6.00]  GHz  [29].  The  WB  search  is  sometimes  rendered  as  WB 
RFSICS  snapshots  by  displaying  a  compilation  of  Wrf  =  36  MHz  sweeps.  Within 
the  RFSICS,  the  collected  signal  is  down-converted  to  an  intermediate  frequency  (IF) 
of  W if  =  70  MHz  and  sampled  at  fs  =  47.5  Msps  using  a  b  =  12-bit  analog-to-digital 
converter  [29].  Finally  the  signal  is  coverted  to  baseband. 

3.3  Simulated  FMCW  Signal 

Generation  of  the  simulated  FMCW  signal  was  based  on  the  model  in  [27]. 
First,  the  noise  was  created  using  an  Additive  White  Gaussian  Noise  (AWGN)  model 
with  a  sample  frequency  fs  =  3.0  GHz.  The  generated  noise  was  bandpass-filtered 
using  a  bandwidth  of  W^pse  =  300  MHz  centered  at  f^mse  =  350  MHz.  Then, 
a  triangular  FMCW  modulated  signal  is  created  with  a  modulation  bandwidth  of 
Wrp  —  300  MHz,  frequency  modulation  period  of  TM=  1.0  /xs,  and  center  frequency 
of  /FM  =  350  MHz.  Finally,  the  noise  signal  was  convolved  with  a  triangular  FMCW 
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Figure  3.2:  Spectral  response  of  the  simulated  FMCW  signal. 

modulated  signal  and  then  was  high  pass  filtered,  removing  the  duplicate  baseband 
signal.  This  produced  a  noise  plus  FMCW  signal,  which  is  the  simulated  FMCW 
signal  with  the  following  parameters  in  Table  3.1.  The  simulated  FMCW  signal 
produced  has  Ns  =  12.0  k  samples,  a  bandwidth  of  Wj^pCW  =  600  MHz  centered  at 
fFMcn  _  7qq  waveform  period  of  7V=  2.0  /xs,  frequency  modulation  period 

of  Tm  =  1.0/zs,  and  is  spectrally  illustrated  in  Figure  3.2. 


Table  3.1:  Simulated  FMCW  Signal  Parameters 


Parameter 

Value 

Sample  Frequency (fs) 

3.0  GHz 

Noise  Power  (Pn) 

-40  dBW 

Bandwidth  (W^cw) 

600  MHz 

Noise  Bandwidth  ( W^pWN ) 

300  MHz 

FM  Bandwidth  ( W 

300  MHz 

FMCW  Center  Frequency 

700  MHz 

Frequency  Modulation  Period  (Tm) 

1.0  /is 

Waveform  Period  (TV) 

2.0  ns 
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Table  3.2:  NTR  System  Components  Configuration 


Hardware 

Number 

Laptop 

S/N:  3MMGRB1 

Transmitter  ID 

AFIT  #4 

A/D  Converter 

AL  83  GTE  8054929 

High  Gain  Antenna 

PCB  Log  Periodic  WA5VJB 

Stand 

n/a 

Cables 

n/a 

3-4  Experimental  NTR  Signal 

The  section  below  describes  the  three  collection  configurations  and  their  collec¬ 
tion  frequencies.  Each  configuration  is  a  slight  modification  of  the  implementation 
monostatic  NTR  system  described  in  Section  2.1  and  is  illustrated  in  Figure  2.2. 
Table  3.2  is  the  NTR  system  hardware  configuration  used  during  the  collections. 
These  configuration  capture  different  levels  of  NTR  signal  “coloration”  from  the 
hardware  with  the  goal  of  extracting  the  NTR  characteristics  and  difference  be¬ 
tween  NTRoff  and  NT  Ron  conditions.  Any  definable  characteristics  in  the  NTR 
Source  configuration,  Configuration  #1,  should  propagate  forward  into  the  other  two 
configurations,  Configurations  #2  and  Configuration  #3.  Also,  any  characteristics 
revealed  in  a  higher  number  configuration  and  not  in  a  lower  number  configuration 
can  be  attributed  to  the  “coloration”  in  that  particular  configuration. 

The  NTR  noise  generator  is  a  WB  thermal  noise  source  [25,30].  For  each 
collection  configuration,  the  corresponding  bandwidth  was  determined  using  WB 
RFSICS  snapshots  that  were  generated  from  the  Wrf  =  36  MHz  NB  filter  scanning 
across  the  fc  e  [0.02,1.6]  GHz  range.  The  resultant  Power  Spectral  Density  (PSD) 
for  collected  NTR  signals  and  each  configuration  are  provided  in  Figure  3.3.  The  re¬ 
sponse  in  Figure  3.3  illustrates  the  NTR  signals  on  a  frequency  and  power  axis  which 
span  /  Rj  [0.02,2.020]  MHz  and  [-110,-40]  dBm.  The  center  frequencies  for  NTR  de¬ 
tectability  analysis  and  demonstration  were  selected  to  coincide  approximately  with 
the  center  of  the  NTR  signal  bandwidth  (Wntr)  and  with  the  lower/upper  roll-off 
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(a)  Config.  #1  NTR  Source  Bandwidth. 


o 

a. 


(c)  Config.  # 3  NTR  Antenna  Bandwidth. 

Figure  3.3:  WB  RFSICS  PSD  snapshot  of  the  collected  NTR  signals 
for  each  configuration. 

regions  of  W^tr,  as  reflected  in  Figure  3.3(a),  Figure  3.3(b),  and  Figure  3.3(c).  The 
selected  fc  for  each  of  the  configurations  are  listed  in  Table  3.3. 


Table  3.3:  RFSICS  Center  Frequencies  (/c)  Used  for  NTR  Collections. 


Configuration  # 

Low  (MHz) 

Mid  (MHz) 

High  (MHz) 

#1 

Direct  NTR  Source  Output 

n/a 

580 

1618 

#2 

Direct  NTR  Tx  Output 

373 

581 

800 

#3 

Over  Air  Tx/Rx 

497 

686 

797 

(b)  Config. #2  NTR  TX  Bandwidth. 


For  collecting  NTR  signals,  there  were  two  conditions  that  were  considered 
for  pair-wise  assessment,  including  NTRqff  and  NTRqn ■  The  NTRqff  condition 
accounted  for  the  collections  including  background- only  noise  responses.  In  this  case, 
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the  NTR  system  was  off  and  not  transmitting.  The  background  noise  consisted  of 
emissions  from  devices  typically  operating  in  an  office  environment,  i.e.,  computers, 
cellular  phones,  power  transformers  in  lights,  etc.  For  Configuration  #1  and  Con¬ 
figuration  #2,  the  background  noise  level  was  low  because  the  RFSICS  was  directly 
connected  to  the  NTR  transmitter.  For  Configuration  #3,  the  collection  receiver  was 
connected  to  a  high  gain  log-periodic  antenna  and  collections  included  environmen¬ 
tal  background  noise.  The  NTRqn  condition  accounted  for  collections  made  while 
the  NTR  system  was  transmitting.  Thus,  the  responses  collected  included  the  NTR 
signal  plus  background  noise  that  was  associated  with  the  collection  configuration. 

3-4-1  Configuration  #1:  Direct  NTR  Noise  Source.  Configuration  #1  is 
illustrated  in  Figure  3.4.  The  RFSICS  was  directly  connected  to  the  NTR  noise 
generator  which  was  thermal  noise  [25,30].  This  configuration  captured  the  charac¬ 
teristics  that  were  truly  inherent  to  the  NTR  noise  generator.  In  this  configuration, 
the  NTR  noise  generator  was  connected  by  a  9.0  ft  cable  to  a  15  dB  attenuator. 
The  attenuator  was  then  directly  connected  to  the  power  divider.  From  the  power 
divider,  the  NTR  signal  was  split  into  the  RFSICS  and  the  ADC,  using  two  separate 
3.0  ft  cables.  The  ADC  then  fed  the  NTR  signal  into  the  spectrum  analyzer  on 
the  laptop,  which  was  used  to  verify  that  the  NTR  system  was  transmitting.  On 
the  other  end,  the  RFSICS  sent  the  NTR  signal  to  another  laptop  that  controlled 
the  RFSICS  and  recorded  the  collections.  For  Configuration  #1,  the  NTR  signal 
Wconfig# i  was  ~  l.T  GHz  and  the  fc  =  for  the  two  36  MHz  NB  collections  frequency 
were  approximately  fc  =  580  MHz  (middle  of  the  H'con/?;g#i),  and  fc  =  1681  MHz 
(high  frequency  roll-off). 

3-4-2  Configuration  # 2 :  Direct  NTR  Tx  Output.  Configuration  #2  is 
illustrated  in  Figure  3.5.  The  RFSICS  was  directly  connected  to  the  NTR  transmitter 
output.  This  configuration  captured  the  NTR  signal  “colored”  by  filtering.  Within 
the  NTR  transmitter,  the  NTR  noise  generator  was  filtered  twice  by  a  low  and  high 
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Figure  3.4:  NTR  Configuration  #1,  direct  connection  from  the  NTR 
Source. 

pass  filter.  The  NTR  signal  was  then  split.  One  end  of  the  NTR  transmitter  output 
was  terminated  by  a  Rl  =  50  O  load  and  the  other  NTR  transmitter  output  was 
connected  by  a  9.0  ft  cable  to  a  15  dB  attenuator.  The  attenuator  was  then  directly 
connected  to  the  power  divider.  From  the  power  divider,  the  NTR  signal  was  split 
into  the  RFSICS  and  the  ADC,  using  two  separate  3.0  ft  cables.  The  ADC  then  fed 
the  NTR  signal  into  the  spectrum  analyzer  on  the  laptop,  which  was  used  to  verify 
that  the  NTR  system  was  transmitting.  On  the  other  end,  the  RFSICS  sent  the  NTR 
signal  to  another  laptop,  which  controlled  the  RFSICS  and  recorded  the  collections. 
For  Configuration  7^2,  the  NTR  signal  WC(mfig# 2  was  ~  400  MHz  [30]  and  the  fc 
for  the  three  36  MHz  NB  collections  frequency  were  approximately  fc  =  581  MHz 
(middle  of  the  Wconfig#2) ,  fc  =  373  MHz  (low  frequency  roll-off),  and  fc  =  800  MHz 
(high  frequency  roll-off). 


NTR  Receiver  Processer 


Figure  3.5:  NTR  Configuration  #2,  direct  connection  from  the  NTR 
Tx  output. 
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3-4-3  Configuration  #3:  Over  Air  Tx/Rx.  Configuration  #3  is  illustrated 
in  Figure  3.6.  The  RFSICS  and  NTR  transmitter  were  connected  to  two  separate 
high  gain  log-periodic  antennas.  This  configuration  captured  the  NTR  signal  “col¬ 
ored”  by  filtering,  the  antennas,  and  air  transmission  background  noise.  Within  the 
NTR  transmitter  block,  the  generated  NTR  noise  signal  was  filtered  twice:  once 
by  a  low  pass  filter  and  once  by  a  high  pass  filter.  After  filtering,  the  NTR  signal 
was  power  divided  with  one  output  terminated  by  an  RL  =  50  D  load.  The  other 
output  was  the  NTR  transmitter  output  connected  by  a  9  ft  cable  to  a  15  dB  at¬ 
tenuator.  The  attenuator  was  then  directly  connected  to  the  power  divider.  From 
the  power  divider,  the  NTR  signal  was  split  into  high  gain  log-periodic  antennas 
and  the  ADC,  using  two  separate  3  ft  cables.  The  ADC  then  fed  the  NTR  signal 
into  the  spectrum  analyzer  on  the  laptop  which  was  used  to  verify  that  the  NTR 
system  was  transmitting.  On  the  other  end,  the  high  gain  log-periodic  antenna  was 
in  close  proximity  (~  2  inches)  to  another  high  gain  log-periodic  antenna  which  was 
connected  to  the  RFSICS  with  a  3  ft  cable.  The  NTR  signal  was  transmitted  and 
received  over  the  air.  The  RFSICS  finally  sent  the  NTR  signal  to  another  laptop, 
which  controlled  the  RFSICS  and  recorded  the  collections.  For  Configuration  #3, 
the  NTR  signal  Wconfig#3  was  ~  400  MHz  [30]  and  the  fc  for  the  three  36  MHz  NB 
collections  frequency  were  approximately  fc  =  686  MHz  (middle  of  the  H/COn/jg#3), 
fc  =  497  MHz  (low  frequency  roll-off),  and  fc  =  797  MHz  (high  frequency  roll-off). 


NTR  Receiver  Processer 


Figure  3.6:  NTR  Configuration  #3,  over  the  air  transmission  and 

receiver. 
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3.5  Post- Collection  Processing 

Figure  3.7  shows  the  post-collection  processing  that  all  collected  data  goes 
through  prior  to  QMFB  processing. 

i 
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J 

Figure  3.7:  Post-collection  processing. 
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3.5. 1  Analysis  Signal  Generation.  For  a  collection  time  period  of  Tc  =  4.0  s 
RFSICS  collection  and  a  resultant  baseband  sample  frequency  of  fs  =  4.75  Msps, 
the  sequence  of  collected  NTR  signal  samples  {R(k)}  contained  Ns  ~  19  million 
samples.  To  make  post-collection  processing  manageable,  a  total  of  Nq  =  2n  «  Ns 
(n  a  positive  integer)  samples  were  taken  from  {R(k)}  to  represent  one  realization 
of  the  NTR  signal.  The  selection  of  the  QMFB  process  input  Nq  =  2n  samples  per 
realization  is  addressed  further  in  Section  2.2. 

3.5.2  Baseband  Butterworth  Filtering.  The  Nq  samples  of  {R(k)}  were 
filtered  using  a  Qth -order  low-pass  baseband  Butterworth  filter  with  a  bandwidth  of 
Wbb  =  18  MHz.  Filtering  the  signal  prior  to  processing  removed  undesired  noise 
components.  The  Wbb  —  18  MHz  was  the  maximum  available  given  an  RFSICS  RF 
bandwidth  of  Wrf  —  36  MHz. 

3.5.3  QMFB  Processing.  The  filtered  Nq  samples  {R(k)}  were  processed 
through  the  QMFB  process  as  described  in  Section  2.2.  This  process  produced 
multiple  non-transformed  QMFB  (Qn)  matrices  for  (n  =  1, 2,  3, ...,  Nl),  where  Nl 
was  the  total  number  of  layers  based  on  (2.6).  A  specific  layer  of  the  non-transformed 
QMFB  Qn  was  chosen  to  be  further  processed. 


3.6  Multiple  Realizations 

Section  4.1  was  based  on  using  one  NTR  realization  with  pre-conditioned  Qctm 
prior  to  graphic  generation.  In  an  attempt  to  increase  visual  features  of  the  NTR 
characteristics,  multiple  realizations  of  the  NTR  signal  were  considered  (Nr  >  1). 
The  Nr  realizations  were  extracted  from  (R(fc)}  and  represented  contiguous  sub¬ 
sequences  of  the  collected  {R(k)}  sequence  containing  Ns  samples.  Therefore,  the 
number  of  available  NTR  realizations  Nr  was  dependent  on  the  length  of  Nr  and  the 
choice  of  Nq.  Figure  3.8  is  an  illustration  of  the  multiple  realization  process  for  the 
pre-conditioned  Qctm  with  additional  pre-conditioning.  The  multiple  realizations 
of  Qctm  were  averaged  to  suppress  any  anomalies  and  to  emphasize  any  details  that 
were  apparent  across  the  realizations.  After  the  averaging,  the  matrix  was  globally 
normalized  to  compare  NTRqff  and  NT  Ron  collection  conditions. 


Figure  3.8:  Detailed  view  of  the  multiple  realizations  of  the  pre¬ 

conditioned  Qctm- 


3.6.1  Pre-Conditioning:  Non-transformed  QMFB.  The  Qn  was  graphically 
represented  using  four  graphic  functions  in  MATLAB®  described  in  Section  2.2  and 
Section  3.7.  These  graphic  representations  were  successfully  used  in  [11,  27]  and 
represent  the  baseline  for  visual  assessment  of  NTR  signals  in  this  research.  However, 
in  an  attempt  to  graphically  extract  more  features,  data  in  the  given  Qn  were  pre¬ 
conditioned  before  graphics  generation. 
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The  following  provides  details  of  various  pre-conditioning  steps  that  were  con¬ 
sidered  in  various  orders  during  the  research.  Each  collection  condition  was  pre¬ 
conditioned  individually,  then  pre-conditioned  pair-wise.  Using  Qn ,  the  first  pre¬ 
conditioning  attempted  was  global  maximum  magnitude  normalization  in  (3.2).  Then, 
the  resultant  energy  was  converted  to  dB  with  (3.1).  Lastly,  pair-wise  thresh¬ 
old/flooring  was  used  for  associated  NTRqff  and  NTRqn  conditions  using  (3.10). 
Multiple  combinations  of  these  three  pre-conditioning  operations  were  attempted  to 
reveal  distinguishable  features. 

3.6.2  Pre-Conditioning:  G-M  Transformed  QMFB.  The  G-M  transforma¬ 
tion  process  provides  an  alternate  means  for  visually  assessing  potential  information 
in  the  non-transformed  QMFB  output,  Qn.  As  adopted  for  this  research  and  depicted 
in  Figure  3.9,  the  G-M  process  has  been  used  for  assessing  the  detectability  of  UWB 
communication  signals  using  a  NB  receiver  [7,8,22,23].  The  G-M  graphics  visual¬ 
ization  in  [7,8,22,23]  was  expanded  under  this  research  using  the  pre-conditioning 
of  Section  3.1  as  illustrated  in  Figure  3.9.  The  need  for  considering  pre-conditions 
was  driven  by  marginal  visual  detectability  using  the  graphic  functions  in  Section  3.7 
with  non-conditioned  G-M  matrices. 

In  the  search  to  obtain  the  “best”  results,  multiple  combinations  of  the  pre¬ 
conditioning  processes  in  Section  3.1  were  attempted  using  trial-and-error.  All  of  the 
considered  combinations  were  qualitatively  assessed  using  the  MATLAB®  graphics 
per  Section  3.7.  The  operations  attempted  are  illustrated  in  Figure  3.10  which  shows 
that  pre-conditioning  was  attempted  both  before  and  after  all  steps  of  the  process. 

3.7  Graphic  Visual  Assessment 

Graphic  visual  assessment  was  performed  for  each  matrix  of  interest,  includ¬ 
ing  non-transformed  QMFB  (Qn),  G-M  transformed  QMFB  ( Q%M '■  Qttmi  Qctmi 
Qssm,  Qc¥m )>  an(l  any  variants  thereof  that  had  have  been  pre-conditioned  per 
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Figure  3.9:  G-M  process. 

Section  3.1.  As  illustrated  in  Figure  3.1,  this  is  accomplished  using  graphic  func¬ 
tions  of  contour,  mesh,  surface,  and  Pcolor  in  MATLAB®.  The  contour  graph  is 
a  2-dimensional  representation  Q(j,  k )  matrix  with  a  time  and  frequency  axis.  The 
mesh  graph  is  a  2-dimensional  representation  Q(j,  k )  matrix  with  an  energy  and  fre¬ 
quency  axis.  The  surf  graph  is  a  3-dimensional  representation  Q(j,  k )  matrix  with  a 
time,  frequency,  and  energy  axis.  The  Pcolor  graph  is  a  2-dimensional  representation 
Q(j,  k )  matrix  with  a  time  and  frequency  axis. 

For  all  cases  considered,  graphics  are  provided  for  pair-wise  assessment  of 
NTRqff  and  NTRqn  conditions.  This  enables  qualitative  assessment  of  “What’s 
different?”  in  the  pair  of  responses.  The  noted  differences  between  background- 
only  and  NTR-plus-background  responses  represent  a  form  of  NTR  signal  detection, 
i.e. ,  the  “presence”  of  differences  suggests  that  the  NTR  signal  has  influenced  the 
composite  response  and  its  “presence”  can  be  declared. 
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Figure  3.10:  Pre-conditioning  operations  attempted  for  G-M  trans¬ 
formed  QMFB  . 


3.8  Wealth  of  Available  Data 

Accounting  for  all  parameters  of  various  processes  depicted  in  Figure  3.1  and 
available  for  analyzing  NTR  detectability,  there  were  a  near-infinite  number  of  sce¬ 
narios  that  could  be  considered  for  analysis.  To  help  quantify  the  number  of  possi¬ 
bilities,  consider  the  following  example: 

1.  Nr  =  1  NTR  realization  (virtually  unlimited  Nr  possibilities) 

2.  Ncc  =  3  collection  configurations 

3-  Nsb  =  9  contiguous,  non-overlapping  NB  sub-bands  required  to  process  entire 
NTR  bandwidth  of  W^tr  ~  400  MHz  [30]  using  the  RFSICS  instantaneous 
bandwidth  of  Wrf  =  36  MHz  (virtually  an  unlimited  number  of  Nsb  as  over¬ 
lapping  sub-bands  are  used) 
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4.  Nq  =  2Nl  =  213  QMFB  input  samples  produces  Nl  =  13  available  Qn  output 
layer  data  and  T-F  matrices  (virtually  unlimited  Nq  =  2Nl  possibilities) 

5.  Nql  —  1  chosen  layer  from  NL  complex  T-F  Qn  (12  of  13  available  Qn  remain 
for  Nq  =  213) 

6-  Nqm  =  4  G-M  transformations  for  Qn 

Accounting  for  parameters  in  #1  though  #6  above,  there  were  Nr  x  Nqc  x  A sb  x 
Nql  x  Nqm  =  108  complex  T-F  matrices  available  for  graphic  analysis.  This  was 
based  on  27  non-transformed  QMFB  Qn  before  G-M  transformation  and  27x4  =  108 
Qttmi  Qctm-i  Qssmi  an(l  Qc¥m  after  transformation.  So  in  total,  there  are  108 
possible  matrices. 

Using  Nql  =  1  chosen  layer,  say  Q 12,  there  were  Nq  =  4  QMFB  T-F  graphics 
and  Nq£ ®  =  2  pre-conditioned  QMFB  T-F  graphic  representations  (12  of  13  available 
Qn  remain).  Applying  G-M  to  given  Q12,  there  were  Nq  l=A  G-M  T-F  graphics  and 
Nq£  m  =  1  pre-conditioned  G-M  T-F  graphic  representations. 

Considering  the  stated  scenario  conditions,  there  were 

Nr  x  Ncc  X  Nsb  X  Nql  X  (Nq  +  NPq£Q  +  N%M  +  NpreGM)  =  297 Graphics  (3.12) 

There  are  297  possible  graphics  for  a  single  realization,  three  collection  config¬ 
uration,  non-overlapping  collection,  and  Nq  =  213  samples.  Considering  multiple  re¬ 
alizations  with  the  same  conditions  stated  above  and  a  collection  duration  of  Tc  =  4s 
and  using  Nq  =  213  samples  results  in  Nr  «  2,  000  realizations  and  more  than  500 
thousand  matrices  and  graphics  would  be  generated.  The  need  for  down-selection 
was  obvious. 

To  complete  the  research  within  available  time  constraints  and  to  ensure  that 
only  useful  graphic  data  are  presented,  results  were  generated  for  a  limited  subset 
of  available  combinations.  From  these  results,  only  a  limited  subset  is  presented 


33 


as  required  to  support  research  conclusions.  The  final  subset  is  based  on  empirical 
down-selection  that  took  place  at  each  step  along  the  research.  The  down-selection 
was  based  on  qualitative  assessment  of  graphics  per  Section  3.7  and  those  repre¬ 
sentations  that  provided  the  most  visually  discernible  features.  The  baseline  for 
down-selection  decisions  is  the  previous  research  of  [27]. 

Summary 

This  chapter  provided  the  details  for  the  detection  process  that  was  built  on 
two  previous  areas  of  research.  [27]  research  provided  the  foundation  and  [7,  8,  22, 
23]  provided  another  avenue  to  view  the  previous  [27]  research.  Detailed  step-by- 
step  results  for  one  realization  of  the  detection  process  and  results  for  the  multiple 
realizations  are  described  in  Chapter  4. 
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IV.  Results 


This  chapter  provides  results  for  NTR  signal  characterization  using  the  qualitative 
visual  assessment  process  outlined  in  Chapter  III.  Section  4.1  presents  the  baseline 
simulation  results  for  the  Quadrature  Mirror  Filter  Bank  (QMFB)  process  using  the 
Frequency  Modulated  Continuous  Wave  (FMCW)  signal  described  in  Section  3.3. 
These  results  are  based  on  QMFB  process  procedures  in  Section  3.5.3  and  are  pre¬ 
sented  for  validation  against  previous  work  in  [27].  This  includes  the  assessment  and 
selection  of  QMFB  process  input  sample  size  (Nq)  and  layer  selection  for  the  NTR 
signal  analysis. 

Section  4.2  and  Section  4.3  provide  a  representative  sequence  of  NTR  results 
obtained  using  the  step-by-step  detection  procedures  outlined  Sections  3.5  and  3.6. 
These  results  were  for  a  single  realization  of  the  NTR  signal  using  a  single  collection 
configuration  as  well  as  fixed  collection  and  post-collection  processing  parameters, 
including:  1)RFSICS  center  frequency  (/c)  and  collection  bandwidth  ( Wrf ),  2)  Nq 
samples  per  NTR  realization,  and  3)  QMFB  layer  number.  This  enabled  the  empiri¬ 
cal  selection  of  the  “best”  graphic  representation(s)  for  visually  distinguishing  signal 
features.  Section  4.2  and  Section  4.3  address  specific  areas  of  interest,  including 
the  effects  of:  1)  signal  collection  and  post-collection  Butterworth  filtering,  2)  using 
the  QMFB  process  (output  designated  as  non-transformed  QMFB,  Qn ),  and  3)  us¬ 
ing  Gronholz-Mims  (G-M)  process  (output  designated  as  G-M  transformed  QMFB, 
Qn)-  Each  section  provides  an  initial  assessment  of  a  single  realization  of  the  NTR 
signal  detectability  using  MATLAB®  graphic  functions  comparing  NTRoff  and 
NT Rqn  conditions. 

Next,  Section  4.4  provides  the  results  using  multiple  realizations  {NR)  of  the 
NTR  signal  using  the  pre-conditioned  Qctm  f°r  each  of  the  three  collection  con¬ 
figurations:  Configuration  #1  in  Section  4.4.1,  Configuration  ^2  in  Section  4.4.2, 
and  Configuration  #3  in  Section  4.4.3.  These  results  were  collectively  assessed  to 


35 


ascertain  the  “best”  graphical  representation (s)  and  parametric  combination (s)  for 
revealing  NTR  signal  features,  i.e.,  those  which  highlighted  characteristic  differences 
and  enabled  discernment  of  NTRqff  and  NTRon  conditions. 

This  chapter  concludes  by  addressing  two  issues  that  represent  opportunities 
for  near-term  research  contribution.  This  includes  preliminary  results  for  assessing 
Spectral  Scalloping  in  Section  4.5.1  and  the  effects  of  Synthesized  Bi-Phase  Coding 
in  Section  4.5.2.  Detailed  analysis  within  each  of  these  areas  could  not  be  completed 
due  to  research  time  constraints. 

4-1  FMCW  Signal  Processing 

The  simulated  FMCW  signal  described  in  Section  3.3  was  analyzed  using  the 
QMFB  process  in  Section  3.5.3.  Note,  the  QMFB  process  requires  an  input  sample 
size  (Nq)  that  is  an  integer  power  of  two,  Nq  =  2Nl,  where  NL  is  an  integer  and  the 
total  number  of  QMFB  layers  for  Qn.  For  the  QMFB  process,  signals  with  a  collected 
sample  size  (Ns)  that  did  not  have  an  integer  power  of  two  have  two  options:  1) 
parse  or  2)  zero-pad  Ns  samples  as  described  in  Section  2.2.  As  adopted  from  [27], 
the  FMCW  signal  generated  has  Ns  =  12,000  samples.  Since  log2(  12,000)  =  13.55, 
Ns  can  either  be  parsed  to  Nq  =  213  samples  or  zero-padded  to  Nq  =  214  samples 
for  processing.  To  focus  on  the  signal,  Nq  =  213  samples  were  first  addressed. 

Using  the  generated  FMCW  signal,  the  QMFB  process  was  performed  to  gen¬ 
erate  multiple  non-transformed  QMFB  (Qn)  matrices  for  each  (n  =  1,  2,  3, ...,  Nl). 
Each  of  these  Qn  matrices  are  graphically  represented  using  the  graphical  functions 
in  MATLAB®  described  in  Section  3.7.  In  particular,  the  four  MATLAB®  graphic 
functions  utilized  are  contour ,  mesh ,  surface ,  and  Pcolor.  These  four  graphical  repre¬ 
sentations  are  illustrated  for  Layer  $4  to  Layer  7^8  of  the  FMCW  signal  in  Figure  2.4 
and  Appendix  A.  The  contour  graphic  was  ultimately  chosen  for  further  analysis 
because  it  visually  presented  the  most  structured  appearance  with  distinguishable 
characteristics  closely  matching  those  of  the  generated  signal.  This  is  illustrated 
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for  the  FMCW  signal  in  Figure  4.1  and  Figure  4.2  which  depict  Qn  as  a  contour 
graphics  for  Nl  =  13  layers.  However,  there  is  no  Q 13  figure  because  Qnl  is  fully 
decimated  to  a  vector.  Consistent  with  previous  work  [11,27],  the  response  in  Fig¬ 
ure  4.3(a),  Layer  7^6,  has  a  triangular  waveform  response  with  a  signal  bandwidth  of 
Wbp  ~  600  MHz,  FM  bandwidth  of  W™  ~  300  MHz,  signal  period  of  Tt  ~  2.0/zs, 
frequency  modulation  period  of  TM  &  l.Ofis,  and  center  frequency  of  fc  ~  700  MHz. 

Figure  4.3  To  properly  choose  the  appropriate  Nq  sample  size  for  further  re¬ 
search,  the  zero-padded  Nq  =  214  samples  of  the  FMCW  were  also  analyzed  using 
the  QMFB  process.  Comparative  Layer  ^6  results  for  Nq  =  2 13  and  Nq  =  214 
samples  are  illustrated  in  Figure  4.3.  When  observing  Figure  4.3,  it  was  evident 
that  the  parsed  FMCW  with  Nq  =  213  samples  (left)  yielded  the  same  signal  infor¬ 
mation  as  the  zero-padded  FMCW  with  Nq  =  214  samples  (right).  As  a  result,  Nl 
is  dependent  on  the  parsed,  signal-only  sample  size  as  in  (2.7). 

The  results  presented  in  this  section  are  consistent  with  the  analysis  provided 
in  [27]  for  QMFB  processing  of  the  FMCW  signal.  Research  in  [27]  concluded  that 
the  contour  graphic  of  Layer  ^6  had  the  same  signal  characteristics  as  presented 
here-validating  that  the  QMFB  process  in  Section  3.5.3  was  properly  implemented 
for  subsequent  application  to  the  NTR  signal.  Thus,  contour  graphics  of  Layer  #6 
based  on  Nq  =  213  samples  per  realization  was  used  for  initial  NTR  investigation. 

4-2  NTR  Signal  Processing 

For  initial  NTR  investigation,  the  collection  scenario  and  post-collection  pa¬ 
rameters  remained  fixed  throughout  the  detection  process  in  Chapter  HI  and  Fig¬ 
ure  3.1.  The  scenario  included  using  collection  Configuration  ^2  with  a  single  real¬ 
ization  of  Nq  =  213  NTR  samples  collected  with  the  RFSICS  using  fc  =  581  MHz 
and  collection  bandwidth  of  Wc  =  36  MHz.  This  section  also  highlights  empirical 
decisions  that  were  made  to  limit  the  number  of  parameters  and  methods  considered. 
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Figure  4.1:  Baseline  FMCW  signal:  Non- Transformed  QMFB  (Qn) 
contour  T-F  Graphics  for  Layer  #1  to  Layer  ^6  using  Nq= 213  samples. 
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Frequency  (Hz)  Frequency  (Hz)  Frequency  (Hz) 


(a)  Layer  #7  (Q7). 
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(c)  Layer  #9  (Q 9). 
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(e)  Layer  #11  (Qn). 


(b)  Layer  #8  (Q8). 


(d)  Layer  #10  (Qio). 


(f)  Layer  #12  (Q 12). 


Figure  4.2:  Baseline  FMCW  signal:  Non- Transformed  QMFB  ( Qn ) 
contour  T-F  graphics  for  Layer  ^7  to  Layer  ^12  using  Ns  =  213  sam¬ 
ples. 
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Figure  4.3:  Baseline  FMCW  signal  with  Ns  =  12,000  samples:  Qq 
contour  T-F  graphics  with  (a)  Nq  =  213  (parsed)  and  (b)  Nq  =  214 
(zero-padded)  samples.  Figure  labels:  A)  fc  =  700  MHz,  B)  W™  = 

300  MHz,  C)  Wbp  =  600  MHz,  D)  Tt  =  2  /rs,  and  E)  Tm  =  1  /is. 

Nq  =  213  samples  was  chosen  because  of  success  achieved  in  [27]  and  re-verified  in 
Section  4.1.  The  two  primary  collection  conditions  are  denoted  as  NTRoff  (only 
background  noise  present)  and  NTRon  (background  plus  NTR  signal  present).  Each 
area  of  interest  in  the  overall  NTR  detection  process  is  detailed  next,  including  the 
effects  of:  1)  signal  collection  and  post-collection  Butterworth  filtering,  2)  using  the 
QMFB  process  (output  designated  as  non-transformed  QMFB,  Qn ),  and  3)  using 
the  Gronholz-Mims  (G-M)  process  (output  designated  as  G-M  transformed  QMFB, 

QGNMY 


4-2.1  Collection  and  Butterworth  Filtering.  Using  NTR  collection  Config¬ 
uration  7^2  in  Section  3.4.2,  received  signals,  R(t),  were  collected  for  the  NTRoff 
and  NTRqn  conditions.  The  collected  signals  were  then  processed  per  Section  3.5. 
As  a  first  step,  the  complex  I&Q  samples  of  R{t)  were  parsed  into  contiguous  se¬ 
quences  (realizations)  of  Nq  =  213  samples.  Figure  4.4  illustrates  temporal  (top)  and 
spectral  (bottom)  responses  that  are  representative  of  one  realization  of  NTRoff 
and  NTRqn  conditions  with  Nq  =  213  samples.  The  amplitude  in  the  NTRqn 
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(a)  Rqff  Realization. 


(b)  Ron  Realization. 


Figure  4.4:  One  realization  for  NTRoff  and  NTRon  conditions 

with  Nq  =  213  samples:  temporal  (top)  and  spectral  (bottom)  re¬ 
sponses. 


condition  is  100  times  greater  than  that  of  the  NT  Rqff  condition.  This  difference 
is  clearly  reflected  in  the  PSD  of  the  responses  as  well. 

For  both  cases,  R(t)  with  Nq  sample  sizes  were  baseband  Butterworth-hltered 
as  per  Section  3.5.2  using  a  bandwidth  of  Wbb  =  18  MHz.  This  produced  the 
Butterworth-hltered  R(t)  responses  and  corresponding  PSD  responses,  illustrated 
in  Figure  4.5.  The  spectral  effect  of  the  Butterworth  filter  is  clearly  evident  by 
comparing  the  higher- frequency  PSD  responses  (regions  above  18  MHz)  in  Figure  4.4 
and  Figure  4.5. 


4-2.2  QMFB  Processing.  Once  filtered  per  Section  4.2.1,  the  resultant 
output  R(t)  with  Nq  =  213  samples  were  QMFB  processed.  This  produced  NL  =  13 
layers  of  non-transformed  QMFB  output  data  that  were  stored  as  complex  I&Q  sam¬ 
ples  in  multiple  Qn  matrices  (n  =  1,  2, ...,  Nl)-  Per  the  discussion  and  justification 
provided  in  Section  4.1,  only  Layer  ^6  results  are  presented  here  for  graphical  genera¬ 
tion  using  the  four  MATLAB®  graphic  functions:  contour,  mesh,  surface,  and  P color 
described  in  Section  3.7.  Figure  4.6  and  Figure  4.7  illustrate  Q$  as  contour  (top-left), 
mesh  (top-right),  surface  (bottom- left),  and  Pcolor  (bottom-right)  for  NTRqff  and 
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(a)  Filtered  Rqff  Realization.  (b)  Filtered  Ron  Realization. 


Figure  4.5:  Post-collection  Butterworth  filtered  NTR  responses  of 

NTRoff  and  NTRon  with  Nq  =  213  samples  temporal  realizations  in 
Figure  4.4. 

NT  Rqn  conditions, respectively.  The  energy  in  NT  Ron  is  1000-times  larger  than 
NTRoff •  The  frequency  ranges  in  the  figures  illustrate  only  the  positive  frequencies 
/  G  [0,18]  MHz  and  the  signal  time  duration  of  Ts  =  0.17  ms.  Initial  assessment 
of  these  graphics  for  all  other  QMFB  layers  revealed  minimal  distinguishable  signal 
characteristics  for  both  NT  Rqff  and  NT  Rqn  conditions.  However,  the  contour 
and  Pcolor  graphics  seemed  the  most  promising  for  subsequent  consideration  and 
were  considered  further. 

As  a  next  step,  Qe  was  pre-conditioned  by  1)  global  maximum  magnitude 
normalizing  the  magnitude  of  Q6  using  (3.1),  2)  converting  the  squared  resultant 
Q 6  energy  responses  to  dB  using  (3.2),  and  3)  pair-wise  flooring  of  the  resulting 
matrix  using  (3.10).  The  pre-conditioned  Qq  was  then  graphically  represented  using 
the  contour  and  Pcolor  graphics  in  Figure  4.8  and  Figure  4.9,  respectively.  This 
ensured  that  resultant  contour  graphics  would  have  the  same  step  size  and  have 
comparable  scales.  Upon  observation,  there  were  only  slight  differences  between 
the  NTRqff  and  NTRqn  responses.  For  the  graphics  considered  thus  far,  there 
was  no  consistent  or  distinguishable  characteristics  in  the  responses  and  minimal 
insight  is  gained  into  detecting  the  NTRqn  condition.  This  motivated  the  need  to 
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(a)  Q 6  Contour  Graphic. 


(b)  Q6  Mesh  Graphic. 


xIO8  xIO  xIO 


(c)  Qq  Surface  Graphic.  (d)  Qe  Pcolor  Graphic. 


Figure  4.6:  Four  graphical  representation  (as  indicated)  of  the  Non- 
Transformed  QMFB  (Q6)  for  Layer  ff6  NTRqff  condition  using  Nq  = 
213  samples. 
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(c)  Qe  Surface  Graphic. 


(d)  Q 6  Pcolor  Graphic. 


Figure  4.7:  Four  graphical  representation  (as  indicated)  of  the  Non- 
Transformed,  QMFB  (Q6)  for  Layer  #6  NTRqn  condition  using  Nq  = 
213  samples. 
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consider  alternatives,  and  the  research  progressed  to  applying  the  G-M  process  after 
the  QMFB  process. 
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Frequency  (Hz) 


(a)  Qe  Response  (dB)  for  NTRqff-  (b)  Qe  Response  (dB)  for  NT  Rqn- 


Figure  4.8:  Pre-conditioned  Non- Transformed  QMFB  Layer  ffQ  (Qe) 
contour  graphics  (dB  scale)  using  Nq  =  213  samples:  (a)  NTRoff  and 
(b)  NT Rqn  Conditions. 


Time  (sec)  x  4  Time  (sec)  x  10J 

(a)  Q6  Response  (dB)  for  NTRoff-  (b)  Qe  Response  (dB)  for  NT Ron ■ 

Figure  4.9:  Pre-conditioned  Non- Transformed  QMFB  Layer  ffQ  (Qe) 
Pcolor  graphics  (dB  scale)  using  Nq  =  213  samples:  (a)  NTRqff  and 

(b)  NT  Rqn  Conditions. 
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4-3  NTR  G-M  Processing 

As  detailed  in  Section  3.6.2,  the  G-M  process  was  used  to  transform  the  out¬ 
put  of  a  given  non-transformed  QMFB  layer  (Qn)  into  four  different  G-M  trans¬ 
formed  outputs.  Given  that  Pcolor  graphics  provided  the  best  clarity  in  previous 
results,  only  Pcolor  graphics  is  considered  henceforth.  Each  G-M  transformation 
produced  complex  output  samples  that  were  stored  in  a  matrix.  These  matrices  are 
identified  here  according  to  the  type  of  G-M  transform  that  has  been  applied  and 
includes:  Qttm  f°r  Temporal- Temporal  Matrix  (TTM)  transformation,  Qctm  f°r 
Cross- Temporal  Matrix  (CTM)  transformation,  Qssm  f°r  Spectral-Spectral  Matrix 
(SSM)  transformation,  and  Qcsm  f°r  Cross-Spectral  Matrix  (GSM)  transformation. 
In  practice,  there  are  actually  two  different  ways  to  generate  correlation-based  Qctm 
and  Qc¥m  results,  one  by  performing  row-wise  cross-correlation  and  another  by  per¬ 
forming  column-wise  cross-correlation  [7,8,22,23]  in  (3.11).  Generally,  the  G-M 
process  produces  a  G-M  transformed  QMFB  represented  as  Q%M ,  and  specific  G-M 
transformations  are  referred  to  as  Qttmi  Qctmi  Qssm-i  ancl  Qc¥m->  where  N  =  layer 
number. 

4-3.1  Single  NTR  Realization.  The  effect  of  the  G-M  process  on  non- 
transfomed  QMFB  output  data  is  illustrated  using  the  same  NTRoff  and  NT  Ron 
realizations  and  Qq  data  from  Section  4.2.2.  Results  for  the  two  non-correlated  Qttm 
and  Qssm  cases  are  presented  in  Figure  4.10,  which  illustrates  Qttm  (top)  and 
Qssm  (bottom)  Pcolor  graphics  for  NTRoff  (left)  and  NT  Ron  (right)  conditions 
as  indicated.  Results  for  two  of  the  four  possible  correlation-based  transformations, 
Qctm  ancl  Qc¥mi  are  presented  in  Figure  4.11  for  row-wise  correlation.  These 
results  illustrate  Qctm  (top)  and  Qc¥m  (bottom)  results  using  Pcolor  graphics  for 
NTRoff  (left)  and  NT  Ron  (right)  conditions  as  indicated.  The  row-wise  results 
are  presented  here  because  they  provided  the  greatest  discernibility.  Note  that  the 
scales  for  NTRqff  and  NTRqn  conditions  in  Figure  4.10  and  Figure  4.11  are  not 
the  same,  and  therefore,  the  results  cannot  be  compared  directly.  However,  among 
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(d)  NT  Ron  Q^m- 


Figure  4.10:  G-M  Transformed  QMFB  Layer  ^6  (Qq)  Pcolor  graph¬ 
ics  using  non-correlated  TTM  (top)  and  SSM  (bottom)  data  with 
Nq  =  213  samples. 


the  hgures  presented  only  the  row-wise  CTM  response  visually  appears  to  have 
a  distinct  pattern.  Therefore,  row-wise  CTM  was  selected  for  generating  Qctm 
and  considered  for  further  evaluation  using  pre-conditioning  to  enhance  the  visually 
discernible  features. 

In  an  attempt  to  enhance  and  extract  features  from  the  Qqtmi  multiple  pre¬ 
configuration  combinations  were  attempted  in  Figure  3.10  from  Section  3.1.  Pre¬ 
sented  below  was  chosen  as  the  “best”  configuration. 


1.  Real  [Qf\  (R) 
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(c)  NTRqn  Qgsm- 


(d)  NTRqn  Qqsm- 


Figure  4.11:  G-M  Transformed  QMFB  Layer  ^6  (Qe)  Pcolor  graph¬ 
ics  using  row-wise  correlated  CTM  (top)  and  CSM  (bottom)  data  with 
Nq  =  213  samples. 
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2.  Row-Wise  Power  Normalization  using  (3.5) 

3.  Calculate  Qttm  from  Section  2.3 

4.  Use  the  Qttm  (C) 

5.  Row-Wise  Cross-Correlation  [7,8]  using  (3.11) 

6.  Take  the  magnitude.  \Qqtm\ 

7.  Threshold/Floor  the  autocorrelation  diagonal  using  (3.9) 

8.  Row-Wise  Maximum  Magnitude  Normalization  using  (3.3) 

Since  cross-correlation  values  are  complex  in  process  ^5  above,  the  absolute  real  and 
magnitude  values  were  evaluated  for  proceeding  processes  #6  through  # 8  above  and 
are  illustrated  in  Figure  4.12.  These  figures  illustrate  the  absolute  real  value  results 
(top)  and  magnitude  value  results  (bottom)  using  Pcolor  for  NTRoff  (left)  and 
NT Rqn  (right)  conditions  as  indicated.  The  magnitude-based  response  seemed  to 
reveal  a  more  distinct  pattern.  Therefore,  the  magnitude-valued  responses  of  the 
preconditioned  Qctm  were  selected  for  further  processing.  However,  discernibility 
between  responses  for  NTRqff  and  NTRqn  conditions  using  the  pre-conditioned 
Qctm  remained  low,  and  additional  improvement  was  sought.  Thus,  the  research 
shifted  from  processing  a  single  NTR  realization  to  evaluating  performance  using 
multiple  NTR  realizations.  This  was  done  to  determine  if  any  “processing  gain” 
could  be  realized  and  provide  increased  visual  discernibility. 

4-4  Multiple  NTR  Realizations:  Pre-Conditioned  Qctm 

Multiple  NTR  realizations  of  the  pre-conditioned  Qctm  were  evaluated  to 
characterize  how  the  NTR  signal  response  changes  over  time.  The  six  sequential 
realizations  of  the  pre-conditioned  Qctm  are  graphically  illustrated  using  Pcolor 
in  Figure  4.13  and  Figure  4.14  for  NTRqff  and  NTRqn  conditions,  respectively. 
These  figures  illustrate  that  the  pre-conditioned  Qctm  response  changes  over  time 
and  that  responses  for  both  collection  conditions  share  some  similar  traits.  To  obtain 
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(c)  NTRoff-  Magnitude  of  Qctm-  (d)  NT  Ron-  Magnitude  of  Qctm- 


Figure  4.12:  Pre-conditioned  Qctm  QMFB  Layer  7^6  (QQ  Pcolor  graphics  using 
the  absolute  real  values  (top)  or  magnitude  values  (bottom)  of  the  CTM  data  with 
Nq  =  213  samples. 
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a  single  comparable  figures,  the  multiple  realizations  of  the  pre-conditioned  Qctm 
were  again  pre-conditioned  by:  1)  calculating  the  multiple  realization  of  Qqtmi 
2)  averaging  the  results  and  3)  global  maximum  magnitude  normalizing  the  data 
using  (3.2).  This  created  multiple  realizations  of  the  pre-conditioned  Qctm  an<i  is 
illustrated  in  Figure  4.15  with  Nq  =  213  samples  and  Nr  =  3000  realizations  for 
both  NTRoff  (left)  and  NTRqn  (right)  conditions.  These  figures  illustrate  some 
slight  differences  between  NTRqff  and  NTRqn  with  a  subtle  pattern  emerging. 
Given  that  these  were  the  most  discernible  features  observed  up  to  this  point  in 
the  research,  multiple  realization  processing  was  carried  into  the  final  phase  of  the 
research,  and  performance  characterized  for  scenario  variation  in  1)  collection  con¬ 
figuration  number,  2)  center  frequency  (/c),  and  3)  number  of  input  samples  to  the 
QMFB  process  (Nq).  The  Nq  values  are  closely  associated  with  the  number  of  re¬ 
alizations  ( Nr )  used  and  the  collected  sample  size  (Ns)  described  in  Section  3.6. 
Therefore,  for  each  Nq,  there  is  a  particular  NR. 

4-4-1  Configuration  #1:  Direct  NTR  Noise  Source.  Results  here  are  for 
Configuration  #1  as  detailed  in  Section  3.4.1  and  depicted  in  Figure  3.4.  The  effect  of 
three  NTR  Nq  sample  sizes  were  analyzed  for  the  two  fc.  Then,  the  “best”  overall 
Nq  and  fc  scenario  were  selected  because  these  values  extracted  the  most  signal 
characteristics  and  highlighted  distinct  differences  between  NTRoff  and  NT  Ron 
conditions  from  Configuration  #1. 

The  multiple  realization  of  the  pre-conditioned  Qctm  f°r  Configuration  #1 
are  graphically-illustrated  as  Pcolor  in  Section  4.4.1. 1,  Section  4. 4. 1.2,  and  Sec¬ 
tion  4.4. 1.3.  These  graphics  were  based  on  using  three  different  Nq  for  fc  =  580  MHz 
and  fc  =  1618  MHz.  When  comparing  the  performance  for  different  Nq  and  fc  val¬ 
ues  and  for  NTRoff  /  NTRon  conditions,  the  differences  appeared  minute,  and 
it  was  very  difficult  to  visually  discern  between  the  two.  However,  one  particular 
Nq  and  fc  scenario  visually  appeared  to  have  a  slightly  more  discernible  difference 
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10  20  30  40  50  60 


(a)  NTRqff  Realization  1. 


(b)  NTRqff  Realization  2. 


(c)  NTRqff  Realization  3. 


(d)  NTRqff  Realization  4. 


Figure  4.13:  Six  sequential  realizations  of  the  pre-conditioned  Qqtm 
P color  graphics  for  NTROFF  using  Ng  =  213  samples. 
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(a)  NT Rqn  Realization  1. 


(b)  NT  Rqn  Realization  2. 


(c)  NT  Rqn  Realization  3. 


(d)  NT  Rqn  Realization  4. 


(e)  NT  Ron  Realization  5. 


(f)  NT  Rqn  Realization  6. 


Figure  4.14:  Six  sequential  realizations  of  the  pre-conditioned  Qc^sm 
P color  graphics  for  NTRon  using  Ns  =  213  samples. 
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(a)  NTRqff'-  pre-conditioned  Qctm-  (b)  NT  Ron'-  pre-conditioned  Qctm- 

Figure  4.15:  Nq  =  213  &  NR  =  3000:  multiple  realizations  of  the 
pre-conditioned  Qctm  f°r  NTRoff  (a)  and  NT  Ron  (b)  conditions 
using  Pcolor  graphical  representation. 

between  NTRoff  and  NT  Ron  conditions.  The  noticeable  differences  were  in  the 
shapes  and  colors,  which  equated  to  energy  levels. 


44.1.1  Configuration  #1:  Nq  =  213  Samples  and  NR  =  3000  Realiza¬ 
tions.  Among  the  two  fc  in  Figure  4.16,  fc  =  1618  MHz  was  chosen  because  for 
NTRoff  in  Figure  4.16(c)  and  NTRon  in  Figure  4.16(d),  these  graphics  had  the 
most  noticeable  differences  between  the  two  collection  conditions  in  the  shapes  and 
colors.  In  NTRqn ,  the  yellow-red  ellipse  on  the  autocorrelation  line  had  slightly 
more  red  color  when  compared  to  NTRoff-  The  middle  yellow-blue-green  vertical 
ellipse  on  the  x-axis  row  number  33  was  darker  green  and  had  more  distinct  yel¬ 
low  areas  with  NT  Ron  which  appeared  to  have  stronger  color  areas  (more  energy) 
compared  to  NTRoff- 

44- 1-®  Configuration  #1:  Nq  =  216  Samples  and  NR  =  600  Realiza¬ 
tions.  Among  the  two  fc  in  Figure  4.17,  fc  =  580  MHz  was  chosen  because  for 
NTRoff  in  Figure  4.17(a)  and  NTRon  in  Figure  4.17(b),  these  graphics  had  the 
most  noticeable  differences  between  the  two  collection  conditions  in  the  shapes  and 
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(a)  NTRqff :  fc  =  580  MHz. 


(b)  NT  Ron-  fc=  580  MHz. 


Figure  4.16:  Configuration  #1,  NTRqff  and  NTRon  multiple  re¬ 
alization  of  the  pre-conditioned  Qqtm  P color  graphical  representation 
with  the  settings  Nq  =  213  samples  and  Nr—  3000  realizations  for  fc 
=  580,  1618  MHz. 
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colors.  The  orange-red  diagonal  ellipses  for  NT  Ron  had  more  yellow  (more  energy) 
connecting  to  the  blue-green-yellow  vertical  ellipse  at  coordinates  [33,20]  and  [33,45]. 


4-4- 1-3  Configuration  #1:  Nq  =  218  Samples  and  NR  =  210  Realiza¬ 
tions.  Among  the  two  fc  in  Figure  4.18,  fc  =  580  MHz  was  chosen  because  for 
NTRoff  hr  Figure  4.18(a)  and  NTRon  in  Figure  4.18(b),  these  graphics  had  the 
most  noticeable  differences  between  the  two  collection  conditions  in  the  shapes  and 
colors.  The  orange-red  diagonal  ellipses  in  NTRon  had  more  yellow  (more  energy) 
connecting  to  the  blue-green-yellow  vertical  ellipse  at  coordinates  [33,20]  and  [33,45]. 
The  center  of  NTRon  had  a  darker  red  area  (more  energy).  Also,  the  orange-red 
diagonal  ellipses  for  NT  Ron  was  straighter  on  the  major  axis  when  compared  with 
NTRqff- 


4-4- 1-4  Configuration  #1  Summary.  Comparing  NTRqff  and 
NTRon  conditions,  there  were  some  subtle  differences  that  seemed  to  occur  in  the 
same  general  response  areas  for  all  Nq  and  fc  considered.  The  NT  Ron  response 
appeared  to  have  stronger  color  areas  and  larger  color  separation  when  compared 
with  NTRoff,  which  equated  to  having  more  energy  in  certain  locations.  Overall, 
Figure  4.18(a)  and  Figure  4.18(b)  for  Nq  =  218  samples  and  fc  =  580  MHz  appeared 
to  provide  the  most  noticeable  differences  between  NTRqff  and  NTRon  conditions. 

4-4-2  Configuration  #2'-  Direct  NTR  Tx  Output.  Results  here  are  for 
Configuration  #2  as  detailed  in  Section  3.4.1  and  depicted  in  Figure  3.5.  The  effect  of 
three  NTR  Nq  sample  sizes  were  analyzed  for  the  three  fc.  Then,  the  “best”  overall 
Nq  and  fc  scenario  were  selected  because  these  values  extracted  the  most  signal 
characteristics  and  highlighted  distinct  differences  between  NTRqff  and  NTRon 
conditions  from  Configuration  #2. 
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(a)  NTRqff :  fc  =  580  MHz. 
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(b)  NT  Ron-  fc=  580  MHz. 
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(c)  NTRqff-  fc  =  1618  MHz. 


^  0.3 

I 

10  20  30  40  50  60 

Row  Number 

(d)  NT  Ron  and  fc  =  1618  MHz 


Figure  4.17:  Configuration  #1,  NTRqff  and  NTRon  multiple  re¬ 
alization  of  the  pre-conditioned  Qqtm  P color  graphical  representation 
with  the  settings  Nq  =  216  samples  and  Nr—  600  realizations  for  fc  = 
580,  1618  MHz. 
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(a)  NTRqff :  fc  =  580  MHz. 


(b)  NT Rqn'-  fc=  580  MHz. 


(d)  NT  Ron  and  fc  =  1618  MHz 
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(c)  NTRqff-  fc  =  1618  MHz. 


Figure  4.18:  Configuration  #1,  NTRqff  and  NTRon  multiple  re¬ 
alization  of  the  pre-conditioned  Qqtm  P color  graphical  representation 
with  the  settings  Nq  =  218  samples  and  Nr=  210  realizations  for  fc  = 
580,  1618  MHz. 
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The  multiple  realization  of  the  pre-conditioned  Qctm  f°r  Configuration  #2 
are  graphically-illustrated  as  Pcolor  in  Section  4.4.2. 1,  Section  4. 4. 2. 2,  and  Sec¬ 
tion  4. 4. 2. 3.  These  graphics  were  based  on  using  three  different  Nq  for  fc  =  373  MHz, 
fc  =  581  MHz,  and  fc  =  800  MHz.  When  comparing  the  performance  for  different 
Nq  and  fc  values  and  for  NTRoff  /  NTRon  conditions,  the  differences  appeared 
minute  and  it  was  very  difficult  to  visually  discern  between  the  three.  However,  one 
particular  Nq  and  fc  scenario  visually  appeared  to  have  a  slightly  more  discernible 
difference  between  NTRoff  and  NTRon  conditions.  The  noticeable  differences 
were  in  the  shapes  and  colors,  which  equated  to  energy  levels. 

4..4..2.1  Configuration  #2:  Nq  =  213  Samples  and  Nr  =  3000  Realiza¬ 
tions.  Among  the  three  fc  in  Figure  4.19,  fc  =  581  MHz  was  chosen  because  for 
NTRoff  hr  Figure  4.19(c)  and  NTRon  hr  Figure  4.19(d),  these  graphics  had  the 
most  noticeable  differences  between  the  two  collection  conditions  in  the  shapes  and 
colors.  In  NTRqn ,  the  yellow-red  ellipse  on  the  autocorrelation  line  had  slightly 
more  red  color  (more  enery)  when  compared  to  NTRqff ■  The  middle  yellow-blue- 
green  vertical  ellipse  on  the  x-axis  row  number  33  was  darker  (less  energy)  in  the 
NTRon ■  The  yellow  diagonal  shape  in  NTRon  had  more  yellow  (more  energy) 
connections  to  the  green  vertical  shape  and  to  the  vertical  walls. 

4. 4. 2. 2  Configuration  # 2 :  Nq  =  216  Samples  and  NR  =  600  Realiza¬ 
tions.  Among  the  three  fc  in  Figure  4.20,  fc  =  373  MHz  was  chosen  because  for 
NTRoff  hr  Figure  4.20(a)  and  NTRon  hr  Figure  4.20(b),  these  graphics  had  the 
most  noticeable  differences  between  the  two  collection  conditions  in  the  shapes  and 
colors.  In  the  middle  light  blue-green-yellow  vertical  ellipse  on  the  x-axis  row  num¬ 
ber  33,  the  green  area  was  more  separated  with  stronger  green-yellow  areas  (more 
energy)  in  the  NTRon  than  NTRoff ■  The  orange-red  diagonal  ellipses  in  NTRon 
had  more  yellow  (more  energy)  connecting  to  the  blue-green-yellow  vertical  ellipse  in 
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(a)  NTRqff :  fc  =  373  MHz. 


(b)  NT  Ron-  fc  =  373  MHz. 


Figure  4.19:  Config.  #2,  NTRoff  &  NT  Ron  multi-realization  of 

the  pre-conditioned  Qctm  Pcolor  graphics  with  the  settings  Nq  =  213 
samples  and  iV^=  3000  realizations. 
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coordinates  [33,20]  and  [33,45].  In  NTRoff,  the  orange-red  diagonal  ellipses  seemed 
to  be  more  straight,  whereas  NT  Ron  had  more  curves. 

4-4-2. 3  Configuration  #2:  Nq  =  218  Samples  and  Nr  =  210  Realiza¬ 
tions.  Among  the  three  fc  in  Figure  4.21,  fc  =  373  MHz  was  chosen  because  for 
NTRoff  hr  Figure  4.21(a)  and  NTRon  hr  Figure  4.21(b),  these  graphics  had  the 
most  noticeable  differences  between  the  two  collection  conditions  in  the  shapes  and 
colors.  The  orange-red  diagonal  ellipses  in  NT  Ron  had  more  yellow  (more  energy) 
connecting  to  the  blue-green- yellow  vertical  ellipse  in  coordinates  [33,20]  and  [33,45]. 
Also,  in  the  three  orange-red  ellipse  on  the  diagonal  axis  for  NTRon ,  there  were 
multiple  darker  red  areas  (more  energy). 

4-4-2. 4  Configuration  #2  Summary.  Comparing  NTRoff  and 
NTRqn  conditions,  there  were  some  subtle  differences  that  seemed  to  occur  in  the 
same  general  response  areas  for  all  Nq  and  fc  considered.  The  NTRqn  response 
appeared  to  have  stronger  color  areas  and  larger  color  separation  when  compared 
with  NTRqff ,  which  equated  to  having  more  energy  in  certain  locations.  Overall, 
Figure  4.21(a)  and  Figure  4.21(b)  for  Nq  =  218  samples  and  fc  =  373  MHz  appeared 
to  provide  the  most  noticeable  differences  between  NTRoff  and  NTRon  conditions. 

4-4-3  Configuration  #3:  Over  Air  Tx/Rx.  Results  here  are  for  Config¬ 
uration  #3  as  detailed  in  Section  3.4.3  and  depicted  in  Figure  3.6.  The  effect  of 
three  NTR  Nq  sample  sizes  were  analyzed  for  the  three  fc.  Then,  the  “best”  overall 
Nq  and  fc  scenario  were  selected  because  these  values  extracted  the  most  signal 
characteristics  and  highlighted  distinct  differences  between  NTRoff  and  NTRon 
conditions  from  Configuration  #3. 

The  multiple  realization  of  the  pre-conditioned  Qqtm  f°r  Configuration  #3 
are  graphically-illustrated  as  Pcolor  in  Section  4.4.3. 1,  Section  4. 4. 3. 2,  and  Sec¬ 
tion  4. 4. 3. 3.  These  graphics  were  based  on  using  three  different  Nq  for  fc  =  497  MHz, 
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(a)  NTRqff :  fc  =  373  MHz. 


(b)  NT  Ron-  fc  =  373  MHz. 


(d)  NT Rqn-  fc  =  581  MHz. 
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(c)  NTRqff-  fc  =  581  MHz. 
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(f)  NT  Ron  and  fc  =  800  MHz 


0.7 

0.6 

0.5 

0.4 

0.3 


(e)  NTRqff-  fc  =  800  MHz. 


Figure  4.20:  Config.  #2,  NTRoff  &  NT  Ron  multi-realization  of 

the  pre-conditioned  Qctm  Pcolor  graphics  with  the  settings  Nq  =  216 
samples  and  iV^=  600  realizations. 
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(a)  NTRqff :  fc  =  373  MHz. 


(b)  NT  Ron-  fc  =  373  MHz. 


(d)  NT Rqn-  fc  =  581  MHz. 


(c)  NTRqff-  fc  =  581  MHz. 
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(f)  NT  Ron-  fc  =  800  MHz 


(e)  NTRqff-  fc  =  800  MHz. 


Figure  4.21:  Config.  #2,  NTRoff  &  NTRon  multi-realization  of 

the  pre-conditioned  Qctm  Pcolor  graphics  with  the  settings  Nq  =  218 
samples  and  NR—  2100  realizations. 
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fc  =  686  MHz,  and  fc  =  797  MHz.  When  comparing  performance  for  different 
Nq  and  fc  values  and  for  NTRoff  /  NTRon  conditions,  the  differences  appeared 
minute  and  it  was  very  difficult  to  visually  discern  between  the  three.  However,  one 
particular  Nq  and  fc  scenario  visually  appeared  to  have  a  slightly  more  discernible 
difference  between  NTRoff  and  NTRon  conditions.  The  noticeable  differences 
were  in  the  shapes  and  colors,  which  equated  to  energy  levels. 

4-4- 3.1  Configuration  # 3 :  Nq  =  2 13  Samples  and  NR  =  3000  Realiza¬ 
tions.  Among  the  three  fc,  fc  =  686  MHz  was  chosen  because  for  NTRqff  hr 
Figure  4.22(c)  and  NTRon  hr  Figure  4.22(d),  these  graphics  had  the  most  notice¬ 
able  differences  between  the  two  collection  conditions  in  the  shapes  and  colors.  In 
NTRqff ,  the  yellow-red  ellipse  on  the  diagonal  line  was  slightly  more  elongated 
and  straight  when  compared  to  NTRon,  where  it  had  some  defining  features.  The 
middle  ycllow-blue-green  vertical  ellipse  on  the  x-axis  row  number  33  had  a  yellow 
line  in  the  NTRqff •  Finally,  the  center,  red  ellipse  was  larger  (more  energy)  with 
a  distinct  shape  for  NTRqn- 

44-3.2  Configuration  #3:  Nq  =  216  Samples  and  NR  =  600  Realiza¬ 
tions.  Among  the  three  fc,  fc  =  686  MHz  was  chosen  because  for  NTRoff  hr 
Figure  4.23(c)  and  NTRon  hr  Figure  4.23(d),  these  graphics  had  the  most  notice¬ 
able  differences  between  the  two  collection  conditions  in  the  shapes  and  colors.  In 
NTRqff,  the  yellow-red  ellipse  on  the  diagonal  line  was  slightly  more  elongated 
and  straight  when  compared  to  NTRon,  where  it  had  some  defining  features.  The 
middle  yellow-blne-green  vertical  ellipse  on  the  x-axis  row  number  33  was  divided 
into  green  areas  in  NTRqn  compared  to  NTRqff,  where  it  was  a  large  green-blue 
area.  The  yellow-red  diagonal  shape  in  NTRon  had  more  yellow  (more  energy) 
connections  to  the  green  vertical  shape  on  the  x-axis  row  number  33.  Finally,  the 
center,  red  ellipse  was  larger  (more  energy)  with  a  distinct  shape  for  NTRqn- 
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(a)  NTRqff :  fc  =  497  MHz. 


(b)  NT  Ron'-  fc  =  497  MHz. 


(c)  NTRqff:  fc  =  686  MHz. 


(d)  NT Ron:  fc  =  686  MHz. 
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(e)  NTRqff:  fc  =  797  MHz. 


(f)  NT  Ron:  fc  =  797  MHz 


Figure  4.22:  Config.  #3,  NTRqff  &  NTRon  multi-realization  of 

the  pre-conditioned  Qqtm  Pcolor  graphics  with  the  settings  Nq  =  213 
samples  and  Nr=  3000  realizations. 
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(b)  NT  Ron-  fc  =  497  MHz. 
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(a)  NTRqff :  fc  =  497  MHz. 
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(c)  NTRqff-  fc  =  686  MHz. 


(d)  NT  Ron-  fc  =  686  MHz. 
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(e)  NTRqff-  fc  =  797  MHz. 
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Figure  4.23:  Config.  #3,  NTRqff  &  NTRon  multi-realization  of 

the  pre-conditioned  Qqtm  Pcolor  graphics  with  the  settings  Nq  =  216 
samples  and  Nr=  600  realizations. 
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4-4-3. 3  Configuration  # 3 :  Nq  =  218  Samples  and  Nr  =  210  Realiza¬ 
tions.  Among  the  three  fc,  fc  =  497  MHz  was  chosen  because  for  NTRoff  hi 
Figure  4.24(a)  and  NT  Ron  hi  Figure  4.24(b),  these  graphics  had  the  most  notice¬ 
able  differences  between  the  two  collection  conditions  in  the  shapes  and  colors.  In 
NTRoff,  the  yellow-red  ellipse  on  the  diagonal  line  was  slightly  more  elongated 
and  straight  when  compared  to  NT  Ron,  where  it  had  some  defining  features.  The 
yellow-red  diagonal  shape  in  NT  Ron  had  more  yellow  (more  energy)  connections  to 
the  green  vertical  shape  on  the  x-axis  row  number  33.  Finally,  the  center,  red  ellipse 
was  larger  (more  energy)  with  a  distinct  shape  for  NT  Ron- 

44.3.4  Configuration  # 3  Summary.  Comparing  NTRoff  and 
NT Rqn  conditions,  there  were  some  subtle  differences  that  seemed  to  occur  in  the 
same  general  response  areas  for  all  Nq  and  fc  considered.  The  NT  Ron  response 
appeared  to  have  stronger  color  areas  and  larger  color  separation  when  compared 
with  NTRqff,  which  equated  to  having  more  energy  in  certain  locations.  Overall, 
Figure  4.24(a)  and  Figure  4.24(b)  for  Nq  =  218  samples  and  fc  =  497  MHz  appeared 
to  provide  the  most  noticeable  differences  between  NTRoff  and  NT  Ron  conditions. 

4-5  Near-Term  Research  Opportunities 

4-5.1  Spectral  Scalloping.  The  NTR  signal  spectrum  in  Figure  3.3  exhibits 
spectral  scalloping  characteristics  similar  to  what  occurs  when  taking  the  FFT  of 
a  truncated  signal  [20].  For  the  NTR  signal  considered  here,  this  is  potentially 
the  result  of  1)  the  NTR  signal  inherently  possessing  such  characteristics,  2)  the 
WB  scanning  of  a  NB  filter  in  the  RFSICS,  or  3)  a  combination  thereof.  For  an 
initial  “Quick-look”  assessment  into  the  cause  of  spectral  scalloping,  the  RFSICS 
was  used  with  NTR  Configuration  7^2  as  a  fixed,  non-scanning  NB  receiver  with 
the  ADC  dynamic  range  fixed  and  the  bandwidth  set  to  Wc  =  5  MHz.  This  NB 
filter  was  sequentially  centered  at  eleven  specific  frequencies  that  corresponded  with 
peak  and  null  responses  in  the  NTR  signal  response  in  Figure  3.3.  The  resultant 
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(f)  NT  Ron-  fc  =  797  MHz 


(e)  NTRqff-  fc  =  797  MHz. 


Figure  4.24:  Config.  #3,  NTRqff  &  NTRon  multi-realization  of 

the  pre-conditioned  Qqtm  Pcolor  graphics  with  the  settings  Nq  =  218 
samples  and  Nr=  210  realizations. 


69 


collected  data  was  averaged,  normalized,  and  its  power  compared  to  the  NTR  signal 
spectrum.  The  results  of  mapping  these  eleven  power  values  to  the  NTR  signal 
response  is  shown  in  Figure  4.25.  As  indicated,  the  resultant  high  and  low  power 
responses  are  relatively  consistent  with  the  NTR  PSD  responses  at  each  of  the  eleven 
frequencies  considered.  This  suggests  that  the  scalloping  inherently  exists  in  the 
NTR  signal.  However,  this  preliminary  finding  requires  additional  investigation  as 
system  calibration  was  not  performed  for  this  initial  assessment.  Given  that  this 
research  focused  on  NTR  assessment  using  non-scanning,  NB  RFSICS  collections, 
the  effects  of  spectral  scalloping  on  the  final  results  in  previous  sections  should  be 
non-existent. 


Figure  4.25:  Top  graphic  is  a  WB  RFSICS  snapshot  of  Configura¬ 
tion  #2  mapped  to  the  bottom  graphic  which  is  the  processed  collec¬ 
tion. 
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4-5.2  Synthesized  Bi- Phase  Coding.  The  effectiveness  of  G-M  process  with 
a  synthesized  bi-phase  coded  NTR  signal  was  assessed  using  the  resulting  detection 
process  and  Qqtm  fr°m  Section  4.2.2.  The  synthesized  bi-phase  coded  NTR  signal 
was  generated  using  the  process  shown  in  Figure  4.26.  To  enable  direct  comparison 
between  uncoded  NTR  and  synthesized  bi-phase  coded  NTR  detectability,  the  col¬ 
lected  that  was  used  to  generate  uncoded  NTRon  results  in  Figure  4.7(a) 

of  Section  4.2.2  was  selected  and  bi-phase  modulated  per  Figure  4.26.  This  particu¬ 
lar  {R(k)}  contained  Ns  =  213  samples.  As  illustrated  in  Figure  4.27,  the  periodic 
bi-phase  coding  waveform  c(t)  was  generated  with  Ng  =  213  samples  and  contained 
four  periods  with  five  “chip”  intervals  per  period.  The  first  symbol  period  in  Fig¬ 
ure  4.27  is  colored  red,  followed  by  three  additional  duplicate  symbols  in  blue.  Given 
that  an  amplitude  value  of  ±1  was  used  for  c(t),  the  bi-phase  coding  did  not  impact 
the  signal’s  energy  or  power  characteristics.  This  can  be  seen  in  Figure  4.28  which 
illustrates  the  PSDs  for  the  uncoded  NTR  signal  (left)  and  the  synthesized  bi-phase 
coded  NTR  signal  (right).  Furthermore,  such  coding  is  generally  undetectable  using 
power-based  or  energy-based  receiver  processing  methods. 


i_ 


Figure  4.26:  Process  for  synthesizing  bi-phase  coded  NTR  signals. 

The  synthesized  bi-phase  coded  NTR  signal  was  processed  and  analyzed  us¬ 
ing  the  final  QMFB  process  presented  in  Section  4.2.2  using  collected  data  from 
Configuration  ^2  with  Nq  =  213  samples  and  Nr  =  1  realization.  The  resulting 
non-transformed  QMFB,  Qq,  contour  graphic  for  the  bi-phase  coded  signal  is  pre¬ 
sented  in  Figure  4.29  and  possesses  virtually  no  distinguishing  features  that  can  be 
associated  with  the  phase-code  characteristics  in  Figure  4.27.  Relative  to  uncoded 
NTRqn  results  presented  previously  in  Figure  4.7(a)  of  Section  4.2.2,  there  are  no 
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1.2  1.4  1.6 


Figure  4.27:  Coding  waveform  c(t)  used  for  generating  the  synthe¬ 

sized  bi-phase  coded  NTR  signal. 


Frequency  (Hz)  x  ^  q? 

(b)  Bi-phase  coded:  NT  Ron 
Figure  4.28:  Synthesized  bi-phase  coded  and  collected  NTR  PSDs. 

discernible  difference.  The  same  Q§  was  then  G-M  transformed  to  produce  Qqtm 
per  Section  4.3.1  and  Pcolor  graphics  were  generated  which  is  illustrated  in  Fig¬ 
ure  4.30(b).  Relative  to  the  corresponding  uncoded  QMFB  result  in  Figure  4.30(a) 
which  is  provided  for  comparison,  the  G-M  process  produces  a  favorable  response  for 
the  phase-coded  NTR  signal,  i.e. ,  there  visually  discernable  islands  of  high  energy 
concentration  along  the  diagonal  that  are  qualitatively  attributable  to  the  phase 
coding  in  Figure  4.27. 

Performance  of  the  G-M  process  with  the  synthesized  bi-phase  coded  NTR  sig¬ 
nal  can  be  viewed  as  either  1)  a  potential  item  of  concern  from  the  NTR  development 
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Figure  4.29:  Contour  Graphic  of  N on-transformed  Qq  synthesized 

bi-phase  coded  NTR  signal  and  NTRon  condition. 
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(a)  Uncoded  NTR  Signal 
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(b)  Bi-Phase  Coded  NTR  Signal 


Figure  4.30:  Pcolor  graphics  for  G-M  transformed  ( Qctm )  QMFB 
Qs  for  uncoded  and  synthesized  bi-phased  coded  NTR  signals. 


perspective  or  2)  a  potential  item  of  promise  from  an  NTR  exploitation  perspective. 
Regardless  of  the  perspective  taken,  the  “presence”  aspect  of  NTR  detectability  has 
been  demonstrated  and  there  is  near-term  opportunity  for  continued  research  in  this 
area. 


Summary 

This  chapter  provided  results  for  the  qualitative  visual  assessment  and  analysis 
of  NTR  detectability.  Various  configurations,  center  frequencies  (/c),  and  number  of 
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Table  4.1:  Summary  of  Multiple  NTR  Realization  Results.  Bold 

entries  denote  “best”  combinations. 


Configuration 

Nq  Samples 

Nr  Realizations 

fc  (MHz) 

2 13 

3000 

1618 

#  1  -  NTR  Source  Output 

2 16 

600 

580 

218 

210 

580 

213 

3000 

581 

# 2  -  NTR  Tx  Output 

2 16 

600 

373 

218 

210 

373 

213 

3000 

686 

#3  -  Over  Air  Tx/Rx 

2 16 

600 

686 

218 

210 

497 

input  samples  to  the  QMFB  process  (Nq)  were  analyzed.  For  collected  unstructured 
NTR  signals,  Table  4.5.2  provides  a  summary  of  the  “best”  fc  and  Nr  scenario 
for  each  configuration  in  boldface  text.  Among  all  the  unstructured  NTR  signals 
collections,  the  best  performance  was  achieved  using  Configuration  ^3  (Over  Air 
Tx/Rx)  with  a  QMFB  input  sample  size  of  Nq  =  218  samples  and  RF  center  fre¬ 
quency  fc  =  497  MHz  because  of  the  fluctuation  of  the  background  noise  and  the 
strong  NTR  signal  produced.  While  this  scenario  provided  the  most  noticeable 
differences  between  NTRoff  and  NT  Ron  conditions  using  P  color  graphics,  the 
overall  discernibility  was  minimal,  and  detection  remained  very  challenging.  For 
synthesized  structured  NTR  signals  created  by  adding  bi-phase  coding  to  a  collected 
NTR  signal  with  identical  configuration  and  collection  parameters  maintained,  there 
were  distinct  differences  between  NTRoff  and  NTRon  conditions.  In  this  case, 
visually-observed  features  were  detectable  and  corresponded  to  the  induced  bi-phase 
code  features.  Ultimately,  the  overall  research  objective  was  met  in  that  the  resul¬ 
tant  G-M  process  was  shown  to  be  useful  for  assessing  and  influencing  both  NTR 
development  (desire  to  maximize  LPD/LPI  potential)  and  NTR  exploitation  (desire 
to  minimize  LPD/LPI  potential)  activities. 
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V.  Conclusion 


This  chapter  summarizes  the  research  activity  and  provides  conclusions  and  recom¬ 
mendations  for  future  research.  This  includes  a  summary  of  the  research  motivation, 
goals,  techniques,  and  results,  which  are  followed  by  a  list  of  recommendations  for 
near-term  follow-on  research  activity.  Supplementary  appendices  and  the  document 
bibliography  are  provided  immediately  thereafter. 

5.1  Research  Summary 

The  concept  of  noise  radar  was  introduced  circa  1950.  However,  its  practical 
implementation  is  only  now  being  realized  in  various  Noise  Technology  Radar  (NTR) 
systems.  The  reduction  in  computational  complexity,  as  brought  about  by  advances 
in  Digital  Signal  Processing  (DSP),  has  generated  near-worldwide  interest  in  NTR 
system  development.  As  development  continues,  NTR  deployment  becomes  more 
inevitable  and  interest  within  the  RF  Intelligence  (RFINT)  community  heightens. 
Thus,  the  RFINT  community  is  becoming  proactive  to  ensure  that  their  arsenal  of 
RF  awareness  tools  are  up-to-date.  The  emergence  of  NTR  systems  will  require  new 
RFINT  capabilities  based  on  non-cooperative,  non-matched  filter  detection  tech¬ 
niques  to  exploit  information  in  emerging  signals  of  interest.  This  research  provides 
an  initial  assessment  of  NTR  signal  detectability  using  such  techniques. 

The  overall  research  objective  was  to  assess  NTR  signal  detectability  using 
qualitative  graphical  visualization.  The  detection  approach  used  a  non-cooperative 
collection  receiver  and  empirical  pair-wise  assessment  of  graphics  obtained  for  NTRqff 
and  NT Rqn  collection  conditions.  The  emergence  of  distinguishable  feature  char¬ 
acteristics  and/or  differences  between  the  two  conditions  dictated  NTR  signal  de¬ 
tectability.  Relative  to  previous  related  research  using  structured  signals  [11,27],  i.e. , 
those  containing  intentional  modulation  features,  the  work  here  addressed  unstruc¬ 
tured  Wideband  (WB)  NTR  signals  that  are  collected  using  a  Narrowband  (NB) 


75 


collection  receiver.  The  signals  of  interest  were  collected  from  a  typical  NTR  node 
in  AFIT’s  Noise  Network  (NoNET)  radar  system.  The  NoNET  system  is  an  inte¬ 
grated  group  of  NTR  nodes  with  centralized  processing  for  data  acquisition,  digital 
processing  and  multi-sensor  image  fusion  [25,30]. 

Two  distinct  NTR  research  veins  are  being  pursued  at  AF1T,  including:  1)  NTR 
developmental  work  aimed  at  improving  NTR  performance  from  a  cooperative  radar 
detection  perspective  and  2)  NTR  exploitation  work  aimed  at  assessing  NTR  ex- 
ploitability  using  non-cooperative  detection  methods-the  focus  of  this  research.  The 
NoNET  NTR  node  signals  are  low  power,  Continuous  Wave  (CW)  signals  that  in¬ 
herently  possess  Low  Probability  of  Detection/Intercept  (LPD/LPI)  characteristics. 
Quadrature  Mirror  Filter  Bank  (QMFB)  processing  was  adopted  here  for  initial 
assessment  of  NTR  signal  detectability.  Consistent  with  earlier  work  that  effec¬ 
tively  used  QMFB  processing  for  detecting  structured  LPD/LPI  waveforms  [11,27], 
detectability  was  determined  through  graphical  visual  assessment  of  energy-based 
graphics  from  non-transformed  QMFB  output  data.  Results  here  for  unstructured 
NTR  signals  revealed  minimal  discernability  between  NTRqff  and  NT  Ron  condi¬ 
tions. 

Improvement  was  sought  using  transformed  QMFB  based  on  previous  AFIT 
research  by  Gronholz  and  Mims  (G-M)  [7,8,22,23].  The  G-M  process  was  effec¬ 
tively  used  for  detecting  LPD/LPI  Ultra  Wideband  (UWB)  communication  signals- 
common  LPD/LPI  characteristics  of  LIWB  and  NTR  signals  suggested  it  was  rea¬ 
sonable  to  consider  G-M  processing  with  NTR  signals.  Relative  to  non-transformed 
QMFB  processing,  results  for  G-M  transformed  QMFB  process  with  NTR  signals 
were  mixed  and  included: 

1.  For  the  collected  unstructured  NTR  signals,  the  best  performance  was  achieved 
using  Configuration  ^3  (Over  Air  Tx/Rx)  with  a  QMFB  input  sample  size  of 
Nq  =  218  samples  and  RF  center  frequency  fc  =  497  MHz  because  of  the  fluc¬ 
tuation  of  the  background  noise  and  the  strong  NTR  signal  produced.  While 
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these  conditions  provided  the  most  noticeable  differences  between  NTRoff 
and  NT  Ron  conditions  using  Pcolor  graphics,  the  overall  discernibility  was 
minimal  and  detection  remained  very  challenging, 

2.  For  synthesized  structured  NTR  signals  created  by  adding  bi-phase  coding  to 
collected  signals  in  1  above,  with  identical  configuration  and  collection  pa¬ 
rameters  maintained,  there  were  distinct  differences  between  NTRoff  and 
NT Rqn  conditions.  In  this  case,  visually  observed  features  were  detectable 
and  corresponded  to  the  induced  bi-phase  code  features. 

The  overall  research  objective  was  met  in  that  the  resultant  G-M  transformed 
QMFB  process  was  shown  to  be  useful  for  assessing  and  influencing  both  NTR 
development  (desire  to  maximize  LPD/LPI  potential)  and  NTR  exploitation  (desire 
to  minimize  LPD/LPI  potential)  activities. 

5.2  Recommendation  for  Future  Research 

This  section  provides  recommendations  for  near-term  research  opportunities 
that  exist  as  a  result  of  this  research.  As  demonstrated  herein,  NTR  signal  de¬ 
tectability  poses  a  significant  technical  challenge  and  there  are  several  areas  that 
could  be  expanded  upon.  Ultimately,  the  goal  of  subsequent  research  would  be  to 
develop  methods  for  improving  NTR  detectability  by  revealing  or  highlighting  more 
substantial  differences  between  NTRoff  and  NT  Ron  conditions.  This  is  particu¬ 
larly  required  when  considering  unstructured  NTR  signal  cases.  Subsequent  research 
activity  could  include: 

1.  Quantitative  Statistical  vs.  Qualitative  Assessment 

While  graphical  visual  assessment  is  useful  when  considering  system  imple¬ 
mentations  having  an  operator  interface,  the  transition  to  automated  detection 
systems  would  generally  require  quantitative  measures  to  characterize  various 
“states”  being  observed  by  the  system,  i.e.,  NTRqff  and  NTRqn  conditions. 
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Such  measures  could  be  provided  through  statistical  characteristic  differences 
that  may  exist  between  the  various  states  considered. 

2.  Alternate  Parameter  Selection  and  Variational  Analysis 

As  noted  in  Section  3.8,  results  and  conclusions  here  are  based  on  a  very 
small  subset  of  the  vast  number  of  available  parametric  combinations.  Addi¬ 
tional  and/or  alternate  combinations  of  the  parameters  considered  here  may 
produce  better  features  that  are  more  readily  apparent  when  graphically  as¬ 
sessed.  Some  of  the  obvious  variations  to  the  research  presented  here  include 
using  a  different:  non-transformed  Qn  layer;  collection  signal  sample  size  (Ns)', 
input  QMFB  sample  size  ( Nq );  collection  receiver  instantaneous  bandwidth 
(Wrf  <  36  MHz  available  using  the  RFSICS);  collection  center  frequencies 
(fc)  that  span  the  full  NTR  bandwidth  (' Wntr )  using  either  contiguous  non- 
overlapping  NB  collections  or  non-contiguous  overlapping  NB  collections. 

3.  Alternate  Collection  Receiver  Architectures 

Alternate  receiver  architectures  include  those  having  a  wider  maximum  in¬ 
stantaneous  bandwidth  than  the  RFSICS  (Wrf  =  36  MHz).  Using  Wrf  = 
100  MHz,  Wrf  =  200  MHz,  ...  ,  would  enable  collection  and  processing  of  the 
full  NTR  bandwidth  (NoNet  Wntr  ~  400  MHz)  using  either  fewer  sequential 
collections  or  a  single  collection  capturing  the  entire  signal  if  Wrf  =  Wntr 
is  available.  In  either  case,  the  instantaneous  collection  of  larger  segments 
of  Wntr  may  be  more  effective  for  capturing  inherent  Time-Frequency  (T-F) 
variation  that  is  not  present  in  shorter  duration,  narrower  bandwidth  collec¬ 
tions. 

4.  Revisiting  Spectral  Scalloping 

As  discussed  in  Section  4.5.1,  the  preliminary  power  assessment  showed  that 
the  observed  spectral  Scalloping  appears  to  be  inherent  in  NTR  Waveform. 
Once  reconfirmed  that  this  is  indeed  the  case,  subsequent  research  could  con¬ 
sider  spectral  scalloping  features  (minimum  null  /  maximum  lobe  spacing,  rel- 
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ative  min-max  levels,  variation  across  multiple  collections,  etc.)  as  potentially 
exploitable  for  characterizing  NTR  signals. 

5.  Modulated  NTR  Signal  Analysis 

NTR  signal  modulation  has  been  used  to  improve  various  aspects  of  NTR 
performance.  There  is  pending  future  research  at  AF1T  that  includes  adding 
some  form  of  modulation  to  NoNet  NTR  nodes  to  improve  Doppler  resolution. 
Thus,  there  may  be  an  opportunity  at  AF1T  to  assess  detectable  structured 
NTR  signals.  Preliminary  results  in  Section  4.5.2  for  the  synthesized  NTR 
signal  with  bi-phase  coding  suggests  that  G-M  transformed  QMFB  processing 
may  be  well-suited  for  this  task. 
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Appendix  A.  Additional  FMCW  Non- Transformed  QMFB  Figures 

Qbi  Qg?  Oi 7  and  Qg 
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Frequency  (Hz) 


(c)  Q 4  Surface  Graphic. 


(d)  Q 4  Pcolor  Graphic. 


Figure  1:  Four  graphical  representation  (as  indicated)  of  the  Non- 

Transformed  QMFB  (Q4)  for  Layer  ff4  NTRqn  condition  using  Nq  = 
213  samples. 
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Figure  2:  Four  graphical  representation  (as  indicated)  of  the  Non- 

Transformed  QMFB  (Q 5)  for  Layer  ^5  NTRqn  condition  using  Nq  = 
213  samples. 
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(a)  Q~i  Contour  Graphic. 


(b)  Q 7  Mesh  Graphic. 


(c)  Surface  Graphic. 


(d)  Qi  Pcolor  Graphic. 


Figure  3:  Four  graphical  representation  (as  indicated)  of  the  Non- 

Transformed  QMFB  (Q7)  for  Layer  ff7  NTRqn  condition  using  Nq  = 
213  samples. 
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(a)  Q 8  Contour  Graphic. 


(b)  Q s  Mesh  Graphic. 


Figure  4:  Four  graphical  representation  (as  indicated)  of  the  Non- 

Transformed  QMFB  (Q8)  for  Layer  # 8  NTRqn  condition  using  Nq  = 
213  samples. 
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Appendix  B.  MATLAB ®  Code 

The  appendix  provides  the  main  MATLAB®  files  used  to  functionally  implement 
the  QMFB  and  G-M  process  as  described  in  Chapter  III  and  used  in  Chapter  IV. 


B.l  QMFB  Process 


Listing  B.l:  Code/QMFB.m 


1  7,  Originally  startpoint.m  by  Dr  .  Pace 
7,  Modified  by : 

7.  Daniel  V.  Atienza  EENG  799  --  Fall  2009 


/  name  =  name  of  orginal  or  parse  file  in  quotes 
6  7«  Tx  =  0(off)  or  1  (on) 


7.  D  layer 
1 


desired  layer  to  be  analyize  for  FreqConstant 
0  for  all  layers 


11  function  Pr oduceLayer s _v2 ( name , Tx , Dlayer ) 
l  clear  all ;  clc ;  close  all  ; 

format  compact 
f ile=name ; 

16  [I  Q  fs]  =  conver t Yt o I Q ( f ile  )  ; 

l _ 

7.  --  Ensures  I  &  Q  are  in  Column  Format  -- 

[v  ,  u]  =  s ize  ( I )  ;  7.  Since  the  QMFBT  work  with  incoming  signals  as 


21 


7.  columns  ,  this  step  is  to  transform  the 

7.  incoming  signal  to  a  column  vector  if  its  a 

7.  row  vector 


if  v==l 


1  =  1';  7.  This  apply  if  the  incoming  signal  is  a  row  vector 
Q  =  Q’;  7.  This  apply  if  the  incoming  signal  is  a  row  vector 


26  end 
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/  --  Formatting  the  Signal  -- 

7. 

T  =  l/fs;  7,  Period  of  the  signal 

31  signal  =  I  +  j  *Q  ;  7«  Forming  the  signal 


°l  disp(’  Length  of 
L=length ( signal ) ; 
7.  disp(’  Number  of 
36  Lc  =  ce il ( log2 (L )  )  ; 
7.  disp(’  Number  of 
z= ( (2 . "  Lc  )  -  L) ; 
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signal  ’  )  ; 

7.  Length  of  the  signal 
Layers  1  )  ; 

7.  It  will  give  the  number  of  layers, 
zeros  in  padding’); 

7.  In  the  case  that  the  signal  length  does  not 
7.  correspond  to  a  power  of  two  the  signal 
7.  will  be  padded  with  zeros  until  the  next 
7.  power  of  2  (2 '(number  of  layers)) 


l  disp ([’ Length  of  signal,  L  =  num2str(L)]) 

7.  disp  ([’  Number  of  Layers,  Lc  =  ’  ,  num2str  (Lc )  ]  ) 

7.  disp  ([’  Number  of  zeros  in  padding,  z  =  ’,  num2str(z)]) 

46 

zz=zeros (1 , z) ;  /  Vector  of  zeros  to  be  padded  to  incoming 

7.  signal  to  do  it  a  power  of  two  in  length 
tt  =  cat(l,  signal,  zz  ’ )  ;  7«  New  signal  vector  to  be  applyed  to 

7.  the  QMFB  resulting  from  concatenate 
51  7.  the  original  signal  with  the  zeros 

NUM_SAMPLES  =  1  engt h  ( 1 1 )  ;  7»  Length  of  the  new  signal 


if  Tx  ==  1 

Tx  =  1  N  ’  ; 

56  else 

Tx  =  ’Bg’  ; 

end 

/ _ 

7.  --  Filtering  the  Signal  with  the  QMFB  --  (qmfb  function) 


61  qmf  b_v2  (Tx  ,  tt  ,  ’  tsinc  ’  ,  f  s  ,  Dlayer  )  ;  7  Here  tlie  function  qmfb  is 

7.  used,  it  will  apply  the  filter  bank  to  the 
7.  signal  with  the  sine  modified  filter 


end 
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7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7 7 7 7 7 7 1 1 7 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7 7 7 7 7 7 7 7 7 7 7 7 7 1 


7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7 7 7 7 7 1 


convertYtoIQ 


7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7 7 7 7 7 7 7 1 1 7 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7 7 
71  7.  Daniel  V.  Atienza  EENG  799  --  Fall  2009 

7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 


7  file  =  converted  data  file  from  RFSICS 


7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 


function  [I  Q  Fsamp]  =  convertYtoIQ ( f ile ) 
76  load (file) 

I  =  real (Y)  ; 

Q  =  imag ( Y)  ; 

Fsamp  =  1/XDelta; 
return 


81 


7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 


7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 


qmf b_ v2 


7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 


86  7  Originally  qmfb.m  by  Dr. Pace 
7  Modified  by : 

7  Daniel  V.  Atienza  EENG  799  --  Fall  2009 


7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 


7  filter  =  tsinc 
91  7  Tx  =  0(off)  or  1  (on) 

7  fs  =  sampling  frequency 

7  Dlayer  =  desired  QMFB  layer  to  anaylized  or  0  for  all 
7  available  layers 
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function  qmfb_v2 (Tx ,f , filter , fs , Dlayer) 


7.  --  Number  of  Layers  to  Compute  -- 

n  =  f  loor  (log2  (length  (f )))  ;  7.  Determine  the  amount  of  layers 

101  7,  from  the  length  of  the  signal 

7.  number  of  layers 
if  Dlayer  ==  0 
N  =  n; 

else 


106  N  =  Dlayer  ; 

end 

save ( ’ N ’ ) ; 

7. _ 

/  -  -  Formatting  the  Signal  --  (to  pass  the  signal  through  the 
111  7.  filter  bank) 

I  (  C 1 : 2  . ~  n]  ,1)  =  f  ([1:2  ,~n])  ; 

out  =  I  ; 

7. _ 


7.  --  Generating  the  Output  Layers  -- 
116  7.  w  =  waitbar  (0)  ; 

7.  Decompose  the  f  unct  ion 
for  lay  =  1:N  7.  layer 

7.  disp(lay)  7«  Show  what  layer  is  been  generated 

flag  =  1;  7.  Flag  used  to  set  up  the  columns  in  the 

121  7.  output  matrix 

7.  waitb ar(l ay /N,w, [’Layer  Progress  - 
7,  num2st r  (  lay  )  ,  ’  /  ’  ,  num2st r  (N)  ]  )  ; 

7.  Reshape  the  output  matrix 
[r  ,  c]  =  size  (  out )  ; 

126  out  =  zeros(r./2,  c.*2); 


for  i 


1:2  .  ~ ( lay  - 1 ) 


7.  column  of  I  (low  to  high) 


[G,H]  =  f  eval  ( f  ilter  ,  I  (  :  ,  i )  )  ;  7»  Evaluate  the  filter 
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*/,  (sine  modified)  over  the 

7.  signal  going  by  columns 
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141 


7,  Setting  the  output  matrix  (Layer)  by  High  Pass  (G)  and 
7.  Low  Pass  (H)  filtering  output 
if  flag 

out  (  :  ,  i.*2-l)  =  H;7»LowPass  Branch  of  QMFB  Tree 

out  (  :  ,  i.*2)  =  G;  7oHighPass  Branch  of  QMFB  Tree 

else 


end  ; 


end  ; 


out  (  :  ,  i.*2-l)  =  G;7oHighPass  Branch  of  QMFB  Tree 

out  (  :  ,  i.*2)  =  H;  7«LowPass  Branch  of  QMFB  Tree 


flag  =  "flag; 


7.  Change  in  the  flag  value 


’/,  Output  Matrix  and  data  to  work  with 
146  I  =  out;  /  Output  matrix  (signal  already  filtered) 

R  =  real(I);  7«  Real  part  of  the  filtered  signal 

Q  =  imag(I)  ;  7o  Imaginary  part  of  the  filtered  signal 

SigSize=  max  (  size  (f  ))  ;  7»Calculates  signal  size 
151  if  Dlayer  ==  0 

7.  Saving  the  data  to  the  same  directory 
save ( [num2str (SigSize)]) ; 

else 

if  lay  ==  N 

156  7.  Saving  the  data  to  the  same  directory 

save ( [num2str (SigSize)]) ; 
end 

end 


161  end ; 

7.  close(w) 


166  tsinc 


1 7. 1. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 7. 1 7. 7. 7. 7. 7. 7.  %  %  7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 7. 1 1 1 1 1 7. 7. 1 1 


7.  Code  by  Dr  .  Pace 


7. 7. 7. 7.  %  7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 


171  function  [dl  ,  cl]  =  tsinc  (cO)  '/,  Sins  Modified  Filter  function 
/  [dl,  cl]  =  tsinc (cO) 

7.  Sine  Modified  Filter  function.  It  use  truncated  sine 
°l  modified  filter  coefficients  to  decompose  the  column 
7,  vector  cO  into  a  (high  frequency)  column  vector  dl  and 
176  7.  (low  frequency)  column  vector  cl. 

7. 

7»  tsinc.su  must  have  been  run  previously  to  create  h .  dat  and 
7«  g.dat  files  that  contain  the  filter  coefficients.  That 
7.  program  determines  how  many  coefficients  are  used. 

181 

l _ 

7.  --  Load  Filter  Coefficients  (Sine  Modified)  -- 

7. 

load  h . dat ;  h  =  h  ’  ; 

186  load  g.dat;  g  =  g’; 

1 _ 

7c  --  Setting  the  Signal  -- 

N  =  length(cO);  7»  Length  of  the  signal  (from  qmfb 

7.  function  "cO"  is  a  column) 

191  pad3  =  length(h)  ./2  +  1;  7.  pad  with  zeros  to  clear  out  filter 

cO  =  [cO;  zeros  (pad3  ,  1)  ]  ;  7.  padding  with  zeros 

7. _ 

7c  Decompose  the  column  vector  cO 

i  =  pad3  :  2  :  (N  +  pad3 -2)  ;  7c  i  will  decimate  by  2 

196 
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7.  compute  cl 
7.  decimate 


54  cl  low  frequency  column  vector 
cl  =  f ilt er ( f 1 iplr (h)  ,  1 , cO )  ; 
cl  =  cl(i); 


201  i  dl  high  frequency  column  vector  dl 
dl  =  f i It er ( f 1 iplr (g)  ,  1 , cO )  ; 
dl  =  dl (i)  ; 

54 _ 


7.  compute  dl 
7.  dec imat e 


B.2  G-M  Process 

Listing  B. 2:  Code/GMcode.m 

54  Willie  H.  Mims  EENG  799  --  Fall  2005/Winter  2006 

54  Updated  and  modified  by: 

/  Daniel  V.  Atienza  EENG  799  --  Fall  2009 


5  '/.function  [CTMRnormR  CTMRabsnormR]  =  .  .  . 

54  GHcode  (x  ,res  ,  Xmatrix  ,  REAL  ,  norm  ,  clr  ) 


54  54 54 54 54 54 54 54 54 54 54 7. 7. 54 54 54 54  54  54  54  54  54  54  54  54  7. 7. 7. 54  54  54  54  54  54  54  54  54  54  54  7. 54 54 7. 54  54  54  54  54  54  54  54  54  54  54  54  54  54  7. 7. 7. 54  54 


7.7.7. 
7.7.7. 
10  7.7.7. 
7.7.7. 
7.7.7. 
7.7.7. 
7.7.7. 
15  7.7.7. 


x  =  R&Q  data  Layer  name 
res  =  IFFT  &  FFT  resolution 
Xmatrix:  1=CTM  or  2=CSM 

REAL:  l=real  or  3=  magnitude  or  anything  else  =  complex 
Norm:  l=Row  power  norm,  2=Col  power  norm 

4=Row  amplitude  norm,  5=Col  amplitude  norm 
anything  else  =  baseline 
clr  =  1  or  0  (min  the  autocorrelation) 


54  54  54  7. 7. 7. 7. 7. 7. 7. 54  7. 54  7. 7. 7. 7. 7. 54  54  54  54  54  54  54  54  54  54  7. 54  54  7. 54  7. 7. 54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54  54 


load (x)  ; 

Nif ft  =  [64 ,  128 ,  256]  ; 

20  Nfft  =  [512,  1024,  2048]; 

Nifft  =  Nifft(res); 

Nfft  =  Nfft(res); 
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55 


M  =  R  +  j  *Q ; 
x  =  M  .  ’  ; 

if  REAL  ==1 

x=  real (x)  ; 
elseif  REAL==3 
x  =  abs (x)  ; 

end 

“/o0/o7o_Row  Normalization  before  processing.  °/» 7«  7  7  °L  °L  °/« 7»  7»  7»  ”/« ”/» 7  7  7  7  7 
if  norm  ==  1 

7.7.77.7.7.7.7.7.  Power  Norm  (Row) 
r  =  x  ; 

rSq  =  r . " 2 ; 

P  =  sum (rSq , 2)  ./size  (r  ,2)  ; 

PinvSq  =  ( 1 . / sqrt (P) )  ; 
rl=  r  ; 

for  t  =  1 : size (r  ,  1) 

rl (t ,  :  )  =  r (t  ,  :  )  . *  PinvSq  (t )  ; 

end 

x  =  rl  ; 

mram 

elseif  norm  ==  2 

7. 7,  mm  7.  Power  Norm  (Col) 
r  =  x  ; 

rSq  =  r . ~ 2 ; 

P  =  sum (rSq  ,  1)  ./size (r  ,  1)  ; 

PinvSq  =  ( 1 . / sqrt ( P ) )  ; 

rl=  r  ; 

for  t  =  l:size(r,2) 

r 1 ( : , t )  =  r ( : , t) . * PinvSq (t) ; 

end 
x  =  rl  ; 
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mmm 


7«  AMP  NORM 

elseif  norm  ==4 

7.7.7._Amp  Normalization  (Row) 

xNorm  =  zeros ( size (x  ,  1)  , size (x  ,  2) )  ; 
rowMax  =  max(abs(x),[],2); 

for  t  =  1 : size (x  ,  1) 

xNorm(t,:)  =  x (t ,:) /rowMax (t)  ; 

end 

x=xNorm ; 
elseif  norm  =  =  5 

7.7.7. _ Amp  Normalization  (Col) 

xNorm  =  zeros (size(x  ,1)  . size (x,2) ) ; 
colMax  =  max(abs(x)  ,  []  ,1)  ; 
for  t  =  l:size(x,2) 

xNorm ( : , t )  =  x(: ,t)/ colMax(t) ; 

end 

x=xNorm ; 

end 


7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 


if  Xmatrix  ==  1;  7.  CTM 
X  =  if f t (x , Nif ft  ,  1)  ; 


type  =  1 CTM  ’  ; 

trans  =  ’ IFFT  1  ; 

elseif  Xmatrix  ==  2  7.  CSM 
X  =  f  f  t  (x  ,  Nf  ft  ,2)  ; 
type  =  ’ CSM  1  ; 
trans  =  ; FFT 1 ; 

end 


7.7.7.  x  =  input  vector /matrix 

7.7.7.  Nfft  =  Nf  ft -point  inverse  DFT  of  vector  x 
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TO  Dim  =  accross  dimention 

7.  NOTE:  M.’  =  (row  ,  col)  =  (freq  ,  time)  therefore  ifft  has  to  be 
1  the  col  to  convert  all  the  freq  to  time 

[Row  Col]  =  size(X); 
clear  CoutC  CoutR 


/o  Column -by-Column  Correlation 


CoutC  =  zeros (Col , Col )  ; 
CoutC  =  X’*X/Row; 

7o  Row -  by -  Row  Correlation 
CoutR  =  zeros (Row , Row) ; 
CoutR  =  X*X 1 / Col ; 

CTMR  =  abs ( real ( CoutR )) ; 
CTMC  =  abs ( real ( CoutC )) ; 
CTMRabs  =  abs (CoutR); 
CTMCabs  =  abs (CoutC); 


7  initialize  matrix 
1  C-by-C  correlation 

7,  initialize  matrix 
‘/,  R-by-R  correlation 


7. 7. 7. 7. 7. 1 7. 7. 7. 7. 7. 7 1 11  Minimize 

if  clr  ==1 


CTMR_min  =  min (min (CTMR) )  ; 

[CTMR_max  CTMR_maxRowLo c ] =max ( CTMR ,  []  ,2)  ; 

CTMRabs _min=min (min (CTMRabs ) ) ; 

[ CTMRabs _max  CTHRabs.maxRowLoc ] =max ( CTMRabs , [] ,2) ; 

for  g  =  1:  length(CTMR) 

CTMR(g , CTMR_maxRowLoc (g) )=CTMR_min ; 

CTMRabs (g, CTMR_maxRowLoc (g) )=CTMRabs_min ; 

end 

end 
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7. 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7.  Normalize 

CTMC_MaxRowVal  =  max ( CTMC ,  []  ,  2)  ; 

CTMCnormR  =  CTMC; 
for  w  =  1 : size ( CTMC  ,  1 ) 

CTMCnormR (w  ,: )  =  CTMC (w  ,  :  )  / CTMC_MaxRowVal (w)  ; 

end 

CTMR_MaxRowVal  =  max ( CTMR  ,  []  ,  2)  ; 

CTMRnormR  =  CTMR; 
for  w  =  1 : size ( CTMR  ,  1 ) 

CTMRnormR (w  ,: )  =  CTMR (w  ,  :  ) / CTMR_MaxRowVal (w)  ; 

end 

CTMCabs _MaxRowVal  =  max(CTMCabs , [] ,2) ; 

CTMCabsnormR  =  CTMCabs ; 
for  w  =  1 : size (CTMCabs  ,  1) 

CTMCabsnormR (w  ,: )  =  CTMCabs (w ,:) /CTMCabs_MaxRowVal  (w)  ; 

end 

CTMRabs_MaxRowVal  =  max(CTMRabs , [] ,2) ; 

CTMRabsnormR  =  CTMRabs ; 
for  w  =  1 : size (CTMRabs  ,  1) 

CTMRabsnormR (w  ,: )  =  CTMRabs (w ,:) /CTMRabs_MaxRowVal (w)  ; 

end 
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